(N 
T— I 

O 
(N 

in 



(N 

in 



Inst ant ons, Twist ors, and Emergent Gravity 



Jonathan J. Heckman^* and Herman Verlinde 

^School of Natural Sciences, Institute for Advanced Study, Princeton, NJ 08540, USA 
^Department of Physics, Princeton University, Princeton, NJ 08544, USA 



Oh' 

Abstract 

(N 
> 

■ Motivated by potential applications to holography on space-times of positive curvature, 

and by the successful twistor description of scattering amplitudes, we propose a new dual 
matrix formulation of A/" = 4 gauge theory on S*^. The matrix model is defined by taking 
the low energy limit of a holomorphic Chern-Simons theory on CP^'^, in the presence of a 
large instanton flux. The theory comes with a choice of S*^ radius £ and a parameter N 
controlling the overall size of the matrices. The fiat space variant of the 4D effective theory 
arises by taking the large scaling limit of the matrix model, with i^i ~ i'^/N held fixed. Its 
. massless spectrum contains both spin one and spin two excitations, which we identify with 

! gluons and gravitons. As shown in the companion paper [1], the matrix model correlation 

5^ \ functions of both these excitations correctly reproduce the corresponding MHV scattering 

amplitudes. We present evidence that the scaling limit defines a gravitational theory with a 
finite Planck length. In particular we find that in the ipi — )■ limit, the matrix model makes 
contact with the CSW rules for amplitudes of pure gauge theory, which are uncontaminated 
by conformal supergravity. We also propose a UV completion for the system by embedding 
the matrix model in the physical superstring. 
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1 Introduction 



While there is httle doubt that string theory provides a consistent theory of quantum 
gravity, it has proven rather difficuh to specify the physical foundations of the theory. Part 
of the issue is that in situations where maximal theoretical control is available, space-time 
is treated as a classical background, rather than as an emergent concept. Related to this, 
the understanding of holography on space-times of positive curvature remains elusive. 

In this paper we propose and develop a new dual matrix formulation of 4D field theory, 
in which the space-time and field theory degrees of freedom simultaneously emerge from 
a large N double scaling limit. The basic idea is to view A/" = 4 gauge theory on 
as an effective low energy description of an underlying topological large gauge theory, 
without any local propagating degrees of freedom. This topological theory takes the form 
of holomorphic Chern-Simons (hCS) theory defined on projective supertwistor space CP^'"^. 
Motivated by the (partial) successes of earlier matrix reformulations of string theory [2H6] , 
and of the twistor string reformulation of gauge theory [7], we consider the hCS theory 
in the presence of a large background flux [S]. We choose this flux to be a homogeneous 
U{N) instanton on S*^ with the maximal possible instanton number k]\f. This special flux 
configuration, known as the U{N) Yang monopolejl] breaks the U{NNc) gauge symmetry 
of the UV theory down to U{Nc), the gauge group of the IR effective gauge theory. 

As we will show, the low energy physics in the presence of the flux is governed by a 
large matrix model, with matrices of size x k^. The guiding principle that deter- 
mines the form and content of the matrix model action is the requirement that the matrix 
equations of motion encapsulate the space of deformations of the Yang monopole, and re- 
duce to the ADHM equations in the appropriate fiat space limit. The finite matrices arise 
because, in the presence of the instanton flux, all charged excitations are forced to move 
in Landau levels. The lowest level is built up from fcjv Planck cells - or fuzzy points - of 
a non-commutative CP^'^, where each fuzzy point represents one elemental instanton. The 
coordinates of the supermanifold become operators, that satisfy a commutator algebra 

[Z^Zt] = 5"^. (1.1) 

The oscillator number is bounded to be less than some integer A^, which translates into an 
upper limit on the allowed angular momentum on the 5*^. 

A striking aspect of this type of non-commutative deformation is that, unlike other more 

^The Yang monopole was actively studied in the context of a proposed 4D realization of the 
Quantum Hall effect in [9]. These authors also found emergent spin 2 modes among the edge 
excitations of the QH droplets, though there are some issues with interpreting these modes as 
connected with an emergent gravitational sector (see e.g. [10] )• Our setup has some similarities 
with the system of [9], but differs in several essential aspects. 
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standard versions of space-time non-commutativity, it preserves all space-time isometries. 
Moreover, because the non-commutative deformation does not directly affect the holomor- 
phic coordinates, the basic geometric link between the twistor variables and 4D physics is 
kept intact. Nevertheless, the instanton density of the background flux introduces a UV 
length scale on the 5*^, which in units of the radius i scales as 



Below this length scale, the physics is described by a holomorphic Chern-Simons theory, 
expanded around a trivial background and without local excitations. At much larger scales, 
however, the collective motion of the instanton background produces rich dynamics that, 
as we will show, takes the form of a quantum field theory with a UV cutoff given by (11. 2p . 
The continuum QFT arises by taking the large N limit while keeping the S"^ radius i fixed. 
One can instead take a combined large and flat space limit, by simultaneously sending 
and i to infinity while keeping ipi fixed. As we will argue, this results in a 4D gravitational 
theory with a finite Planck length. 

The emergence of gravity from an a priori non-gravitational dual system is a remarkable 
phenomenon that deserves further exploration. At a superficial level, our matrix model 
looks like a regulated version of the twistor string theory introduced by Witten as a dual 
description of A/" = 4 gauge theory. Twistor string theory also gives rise to gravitational 
degrees of freedom in the form of conformal supergravity. However, our set-up differs in 
several essential aspects from the twistor string. Most importantly, our theory comes with a 
natural length scale that breaks conformal invariance. Moreover, the emergence of gravity 
has a clear geometric origin via the non-commutativity of twistor space. 

The non-commutativity (11. ip of the twistor coordinates causes gauge transformations 
to act like diffeomorphisms. Correspondingly, a non-zero gauge field background on non- 
commutative twistor space translates into a deformation of the emergent 4D space-time ge- 
ometry. This relationship can be made precise, via a direct variant of Penrose's non-linear 
graviton construction of self-dual solutions to the 4D Einstein equations. In a companion 
paper |1] we use this insight to evaluate correlation functions of bilinear matrix observ- 
ables, and show that quite remarkably, the end result reproduces MHV graviton scattering 
amplitudes of ordinary 4D Einstein (super)gravity. 

The matrix model is the low energy approximation to a more complete UV theory. In 
the first part of the paper, the UV theory is taken to be the hCS theory, or topological B- 
model open string theory, defined on super twistor space. In the last section, we propose an 
alternative UV completion in terms of type IIA superstring theory. The brane configuration 
that we consider is given by a DO-brane bound to a stack of D8-branes and an 08-plane. 
The background geometry is taken to be 10 non-compact space-time dimensions. In the 
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Figure 1: The twistor matrix model captures the low energy physics of holomorphic Chern- 
Simons theory on CP^I^ in the presence of an instanton background. The ingredients of 
the matrix model have a natural connection with the ADHM construction of instantons. 
Correlation functions of the matrix model compute amplitudes of a 4D space-time theory, 
which in the fiat space limit coincide with scattering amplitudes of gluons and gravitons. 

presence of a background magnetic flux, the low temperature limit of a gas of fermionic 
0-8 strings reproduces the degrees of freedom and action of the twistor matrix model. An 
intriguing aspect of this construction is that it (potentially) gives rise to a 4D theory coupled 
to gravity, without first having to choose some Kaluza-Klein compactification manifold. 

The rest of this paper is organized as follows. In section [2] we review some basics 
of twistor geometry. Next, in section [3] we introduce holomorphic Chern-Simons theory 
with the Yang monopole flux and explain how it naturally gives rise to a non- commutative 
deformation of twistor space. This is followed in section lUby the presentation of the twistor 
matrix model, and the gauge currents of the defect system. Section |5]is devoted to matching 
a particular limiting case of the matrix model to the ADHM construction of instantons. We 
study the continuum limit of the matrix model in section [61 Section [7] discusses the sense in 
which the theory describes a quantum theory of gravity, and in section |8] we provide a UV 
completion of the matrix model based on thermodynamics of a stringy brane construction. 
We present our conclusions and avenues of future investigation in section O 

2 Twistor Preliminaries 

In this section we review some preliminary aspects of twistor geometry, as well as the 
connection between hCS and self-dual Yang Mills. We also discuss quantization of twistor 
space. 
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2.1 Twistor Correspondence 



Twistors were introduced by Penrose as a natural geometric description of massless field 
theories [II1[12] (see e.g. [T3l[Tl] for reviews). The construction starts with the introduction 
of spinor helicity variables. In complexified Minkowski space, light-like momenta Paa can be 
represented as a product paa = TCaT^a of a chiral and an anti-chiral spinor. Next, one notices 
that given the spinor tTq and a space-time point x"" on complexified Minkowski space, one 
can define a corresponding two component complex spinor cu" via 

UJ^ = ix^'^TTa. (2.1) 

This relation is invariant under simultaneous complex rescaling of the two spinors TTa and w". 
One thus obtains a set of homogeneous coordinates Z"" = (w", tt^) on projective twistor space 
PT*. Similarly, one can introduce a dual twistor space, with coordinates Zjj = (Tia,u°'), 
denoted by PT,. Geometrically, PT, and its dual both define a three-dimensional complex 
projective space CP^. We will raise and lower two component spinor indices with the help of 
the corresponding two component epsilon symbolsjfl The associated antisymmetric spinor 
inner products are denoted via 

(ttiTTs) = e^ViaVTaft, [^1^2] = £ab^i002 ■ (2-2) 

Equation (12. ip combines two complex linear relations in CP'^. For a given point x in 
complexified space-time, it defines a CP^ , which we will call the twistor line associated 
with X. The space x CP^ with coordinates (x, vr) is known as the "correspondence space" . 
This admits a map to CP^ via = ix°'°"Ka- 

A convenient presentation of complexified, conformally compactified Minkowski space 
is given as the zero locus of the Klein quadric in 



e^p^iX'^^X^' = . (2.3) 

Here the X"'^ = —X^"" are the six homogeneous coordinates of CP^. The constraint (12. 3p 
is automatically solved by introducing a pair of points in twistor space, with homogeneous 
coordinates Z"" and , via 

_ (2.4) 

These homogeneous coordinates satisfy: X^"'^ Z'^^ = 0. For a given X°'^ these equations are 
solved by all points on the line through Z and W. 

Since the X"'^ are homogeneous coordinates, they are only sensitive to the conformal 



^Our conventions are as follows: cj^ = uj"'e^j^, vr'' = vTae" 
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structure of space-time. Conformal symmetry is broken by designating a choice of two index 
anti-symmetric bi-twistor, called the infinity twistor, which we denote by Iai3 = —Ij3a- The 
'inverse' bi-twistor is denoted by 1°"^ = ^e°'^'^^I^s- The infinity twistor allows us to raise 
and lower the index of the twistor coordinates Z° via = laf^Z^^. Furthermore, it defines 
an anti-symmetric pairing (ZW) between two twistors Z and W via 

{ZW) = lafsZ'^W^ (2.5) 

As the name suggests, the infinity twistor selects the twistor lines that map to space-time 
points at infinity: a pair ZW describes an infinite point if (ZW) = 0. Deleting the locus 
(ZW) = 0, we can introduce affine coordinates x"^, and their duals Xafs, viM 

Z^^'W'^^ . 1 g 

= (^ZW) ' ^"'^ ^ 2^ai3-fSX^ (2.6) 

Raising and lowering of a single index is accomplished via the infinity bitwistor. These co- 
ordinates satisfy the relations x"-^x.y^ = —x^^, and Xa-yX"'i3 = —Xa/3. The dual coordinates 
Xai3 act as projection matrices onto the twistor line associated to the space-time point X: 

Xo^^U^ = (2.7) 

which is solved by all points U = aZ + bW on the twistor line through Z and W. 

The form of the infinity twistor depends on the choice of space-time signature and 
metric. In most of this paper, we will work in Euclidean signature. Euclidean twistor space 
can be obtained from CP^ by introducing a map a which acts on the coordinates as: 



{Z\Z^,Z',Z^) A (z\-Z\z\-Z'^ (2.8) 



where Z denotes the complex conjugate of Z". The point Z" and its image W°' = a {Z°') 
define a corresponding space-time point via equation (12. 4p . None of the resulting twistor 
lines intersect. This is why CP^ can be described as an 5^ fibered over S^. 

The infinity twistor of the euclidean four sphere of radius i can be chosen of the form 

(2.9) 




■^Note that pairs ZW with {ZW) = indeed map to points in complexified space-time with 
j-a0 — S^, this infinity locus does not intersect with the real four sphere. In the flat space 

Hmit, ^ — )• 00, the infinity twistor degenerates, and selects VTa = as the twistor line associated 
with space-like infinity. 
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When 7 = 1, we have a round 5"^ with isometry group 5*0(5). When 7 = 0, we have the 
infinity twistor for Minkowski space. Decomposing the twistor into two component spinors 



Z" = (a;^v^J, Z^ = -(^^;), (2.10) 



the twistor inner product takes the form 



= ^(vTiTTs) + j[uJlUJ2] , (2.11) 

with (..) and [..] the spinor pairings (12. 2p . Note that, unhke for Minkowski space, the 
infinity twistor (12. 9 p for non-zero 7 is invertible. 

The correspondence between the two twistor fine equations (12. ip and (12. 7p is made as 
follows. The dual space-time coordinates can be parameterized with the help of five 
coordinates constrained to live on a round 5^ of unit radius, via 




y^yA = i. (2.12) 



^(1 -?/5)edfe. 



Here = ^y^iic^^)^ with = (a*, — ^1) the usual Pauli matrices. With this parametriza- 
tion, the four component twistor line equation (12. 7p reduces to the standard two-component 
twistor line equation (12. ip with 

-b = ^- (2.13) 

The fiat space limit amounts to zooming in on the south pole region of the near ^ —1. 
In this limit, the remaining four S"^ coordinates y'^ become identified with the fiat space 
coordinates x'^. 

The geometry of twistor space also extends to theories with supersymmetry. Here 
we focus on theories with A/" = 4 supersymmetry. This is achieved by supplementing 
the four bosonic coordinates by four fermionic coordinates i = 1,...,4. The resulting 
supertwistor space CP^'^ is a Calabi-Yau supermanifold. The homogeneous coordinates of 
CP^I^ will be denoted as = {Z'^ltp'). Chiral Minkowski superspace is then described by 
coordinates [x"^,9^^). One can also introduce coordinates on Minkowski superspace given 
as (x""^, x"", 6*" j, 6'!^?') subject to the constraint equations: 

^aa ^ xi^<^±eifi^ (2.14) 

so that the independent coordinates on complexified superspacetime M^'^^ are (x, 6'" j, 6'!f ). 
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Chiral and anti-chiral superspace are then parameterized by {x"f',6'^^) and {x'!^,6'^^), re- 
spectively. Let us note that in the complexified setting, 9+ and ^_ constitute independent 
coordinates. A bosonic line in supertwistor space is specified by: 

UJ^ = ixt^7T^ ■ ^^ = e'^TTa (2.15) 

which are invariant under re-scalings of the superspace coordinate , and so define ho- 
mogeneous coordinates on projective supertwistor space CP^'^. Equations fl2.15p define a 
bosonic CP^''^ in CP^'^. The correspondence space is now parameterized by coordinates 
(x,e_,7r) for C^is X CP^I°, with the map to CP^'^ specified by equations f l2J[^ . 

There is a natural extension of the infinity twistor to a super bi-twistor, X/j. In the 
case of S"^'^, the super bi-twistor is: 



ij 




where rjij is a symmetric tensor. 

2.2 Holomorphic Chern-Simons 



\ 

e"'' (2.16) 
Vij J 



Two of the most notable successes of twistor theory are the construction of anti-selfdual 
Yang-Mills equations, and the application of twistor techniques to the study of perturba- 
tive scattering amplitudes in A/" = 4 gauge theory. Central to both applications is the 
link between 4D (self-dual) gauge theory and holomorphic Chern-Simons (hCS) theory on 
twistor space. Following Witten, we will consider the hCS theory as the field theory of the 
physical modes of the B-model topological open string theory, attached to suitable twistor 
space-filling D-branes. 

Supertwistor space CP'^'^ is a Calabi-Yau supermanifold, making it possible to introduce 
a topological B-model. The corresponding open string subsector defines a decoupled theory. 
The physical degree of freedom is a (0, 1) component ^ of a gauge connection on CP'^''^. 
The action is given by the holomorphic Chern-Simons form 

Shcs = I n AtifAdA+'^AAAAA) (2.17) 
The holomorphic (3|4) form Q for the Calabi-Yau space CP^'^ is given by: 

n = ^^z^^dz^dz^dz^ ^JihLd^p'dip'diij^diijK (2.18) 



8 



where the and tp'^ for a = 1 , . . , 4 and i = 1 , . . . , 4 are respectively bosonic and fermionic 
homogeneous coordinates for CP'^'^. The measure is invariant under PS'L(4|4)Q 

The critical points of Shcs correspond to gauge bundles with vanishing (0, 2) curvature: 

J^^'^'^^ = dA + AaA = 0. (2.19) 

In other words, they are holomorphic bundles whose complex structure is compatible with 
the complex structure of the ambient CP^'"^. Via the Ward correspondence [15] (see also 
[T6l[T7]), there is a special class of such holomorphic bundles which describe anti-self-dual 
instantons on S*^. The basic equivalence is between the following two types of mathematical 
objects: 

• Anti-self dual Yang-Mills connections of GL{n,C), defined on the four sphere S"^ 

• Holomorphic rank n vector bundles on CP'^, that are trivial on every twistor line 

The instanton number of the gauge bundle on S*^ corresponds to the second Chern class 
of the vector bundle on CP^'^. Ward's construction represents the analog for gauge fields 
of Penrose's 'non-linear graviton' construction of the anti- selfdual Einstein equations [TS] . 
For more details of the Ward correspondence, we refer to the standard literature [T5|[T6]. 

The connection between the variables of twistor space, and super Yang-Mills theory can 
also be seen through the expansion of A in the fermionic coordinates: 

A = a + + jf^^<P^, + ^^V^'V^VNfo-fc] + j^^'^'^^'^Sijki] (2.20) 

The variable ip has homogeneity -|-1. From this we conclude that a is a degree zero (0, 1) 
form on bosonic twistor space, and 6 is a degree —4 (0, 1) form. Via the correspondence 
between homogeneity and space-time helicity (reviewed in footnote 11), these components 
of the expansion are to be identified with the components of the A/" = 4 vector multiplet. 

Expanding out the components and performing the fermionic integrations, the holo- 
morphic Chern-Simons action reduces to the the following action defined over the bosonic 
twistor space CP^: 

•S'hcs = / Vl' Miih A (da + a A a) + (d + a)(l)ij) + fermionic terms. (2.21) 

The bosonic holomorphic three-form Q' is an element of Q^^'^\0 (4)), while the integrand is 
an element of il^^'^\0 (—4)). The equation of motion for b selects a choice of a such that the 

"^To be slightly more precise, the dip entering here is more appropriately viewed as a measure 
on the Berezinian in the fermionic directions (?]• 
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(0, 2) component of the corresponding field strength vanishes. This defines a holomorphic 
bundle over CP^, which includes all anti-self-dual instantons on S^. Hence the hCS theory 
on CP'^'^ is physically equivalent to the anti-self-dual sector of A/" = 4 gauge theory [7], 
provided that A is restricted to be flat over each twistor line. 

In the influential paper [7], Witten proposed to extend this relationship to the com- 
plete A/" = 4 gauge theory. The construction, motivated by Nair's observation that MHV 
amplitudes are localized on twistor lines, involved supplementing the hCS theory with 
Dl-instantons, wrapping holomorphic curves inside CP^'^. In this paper, we will follow a 
different procedure: we will consider the pure holomorphic Chern-Simons theory, without 
adding any additional degrees of freedom. Rather than introducing Dl-instantons by hand, 
we will study the perturbative dynamics around a non-trivial instanton background for the 
gauge field A on CP^'^, with a very large instanton number. The effect of this flux back- 
ground is quite similar to adding Dl-brane^, in that it leads to additional defect degrees of 
freedom that localize on the twistor lines. But there are important differences with Witten's 
proposal. Most notably, our setup introduces a UV scale, via the size of the constituent 
instantons that comprise the total flux background. 

Two additional comments are in order. First, note that not all holomorphic vector 
bundles on CP^'^ correspond to instantons on S*^: the critical points of the hCS theory 
also include vector bundles which do not trivialize over twistor lines. Mathematically, such 
bundles leads to additional data beyond that specified in the usual ADHM construction 
of instantons on S*^]^ Secondly, the reduction from a string field theory to the topological 
B-model may also allow for certain deformations that can not be captured purely in terms 
of purely holomorphic data. In particular, like in type II theory, one can choose to turn on 
a 2-form S-field background, which deforms the string worldsheet action. The 2-form field 
acts like a U{1) magnetic field, under which the two open string end points are oppositely 
charged. Upon taking the decoupling limit, each end point is forced into lowest Landau 
level orbits. From the point of view of the hCS theory, turning on this 2-form i?-field 
implements a non-commutative deformation of twistor space. 

2.3 Twistor Quantization 

Twistor space and its dual space naturally combine into a quantized phase space |12j . 
Indeed, the Heisenberg commutation relation between positions and momenta implies that 

^Indeed, lower dimensional branes embedded inside of higher dimensional branes dissolve into 
topologically non-trivial flux configurations of the higher dimensional gauge field. 
^See for example [19] for some discussion of this extension in the math literature. 
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the twistors = (cu", tIq) and dual twistors Zp = (ttj,, co^) satisfy the commutation algebra: 

[Z'^,Z,] = m^ ^ (2.22) 

The Hilbert space realization of this commutator algebra provides a representation of 
s/(4, C), the algebra of the complexified conformal group in four dimensions, generated 
by respectively the translations, conformal boosts, Lorentz rotations and dilatation 



Pda — T^aT^a , Kad — ^a^a i 



Jdb = ''^(a^h) ' Jab = ^(aT^h) 



D=]^{naUJ^ -u'^TTa). (2.23) 



The passage to a real space-time signature is achieved by a convention for hermitian 
conjugation of the variables and the corresponding Zp variable, or equivalently, a choice 
of realization of s/(4, C). The three signatures (+ H ), ( — h ++) and (+ + ++) corre- 
spond to the realizations of s/(4, C) respectively by s/(4, M), sm(2, 2) and sm(4). Hermitian 
conjugation in the first two cases acts via 

Tia and G M (2.24) 
\ Biid he iM (2.25) 

T^a J 

The reality condition on fi is enforced by considering how hermitian conjugation acts on 
the canonical commutator relationslll For euclidean signature, hermitian conjugation relates 
the TT and uj spinors via 

sm(4) : ul = na, n^"- = and ^ G M. (2.26) 

The resulting commutation relations are: 

[uj\ujl]=n5l ■ [TT^vrt^] =M^ (2.27) 

Note that daggering an oscillator changes an upper index to a lower one and visa versa. 
Hermitian conjugation fl2.26p acts on the conformal symmetry generators as: 

''The above reality conventions may be deformed by including relative complex phases. This 
provides a means to analytically continue from one space-time signature to another. 



sm(2,2) : 
s/(4,M) : 



^1 = ^a, T^l 
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Figure 2: Depiction of the flat space limit, corresponding to a stereographic projection of 
the 5*^ to M^. In the oscillator algebra this corresponds to a Wigner-Inonii contraction of 
so(5) to the Euclidean Poincare algebra. 

The Euclidean theory is therefore most naturally viewed as a radially quantized theory on 
an S*^. Working with respect to this realization, note that in su{A) ~ so(6) the infinity 
bitwistor (12.91) transforms as a vector under so(6). The generators which leave it invariant 
span a usp{4) ~ so(5) subalgebra, which we identify with the symmetries of the four sphere. 
The Hermitian so(5) generators corresponding to motion along a great circle are: 

Vaa = Paa + i^'"- (2.29) 

Starting from S*^, we can take the fiat space limit by sending the radius of the sphere to 
infinity and zooming in on a local fiat region of the sphere. This procedure is equivalent to 
performing a Wigner-Inonii contraction, and amounts to rescaling the dual oscillators as: 

^1 = 7"^7ra, = l^a (2.30) 

where the parameter 7, the same one that appears in the 5*^ bitwistor (12. 9p . is being sent to 
zero. In this limit, the so(4) subalgebra spanned by J and J stays fixed, while the Hermitian 
translation generators become P + '^'^K. The contraction to the Euclidean Poincare algebra 
then proceeds via 7 — )■ 0. See figure [2] for a depiction of the fiat space limit. 

The extension to space-time supersjTumetric theories is achieved by including four 
fermionic oscillators with j?/'*, ■i/'j} = which then provides a representation of 
ps/(4|4). In what follows we shall often leave implicit the extension to the fermionic case. 
The supersymmetric hermitian charges are the symmetries of the supersphere S"^'^. 

Finally, we emphasize that the above discussion follows the standard viewpoint of twistor 
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quantization [12]. In particular, up to this point, the introduction of quantized twistor 
space does not lead to any a priori space-time non-commutativity. A key feature of twistor 
theory is that space time is a derived notion, obtained via the correspondence at the level 
of holomorphic geometry. Since (for any of the three possible space-time signatures) the 
holomorphic twistor coordinates = (w", tTo) all commute among each other, the notions 
of twistor lines and space-time points are essentially unaffected by the quantum relation 
with the dual twistors. Hence in most applications of twistor methods - for example, the 
construction of instantons and the study of scattering amplitudes - one does not need take 
into account the quantum nature of the twistor coordinates. 

3 Instantons and Fuzzy Twistors 

In the following sections, we will introduce a twistor matrix model, which in a large 
scaling limit reproduces the 4D continuum theory. The matrix model is to be viewed as 
the low energy limit of a more UV complete theory given by ordinary holomorphic Chern- 
Simons theory on ordinary supertwistor space. By expanding around a flux background 
which retains the symmetries of the 4D space-time, we obtain a low energy effective theory. 
Motivated, in part, by the correspondence with the 4D Quantum Hall effect studied in [H], 
the gauge field that we will turn on is the twistor lift of the Yang monopole [20l[2T] 




The Yang monopole is a non-abelian generalization of the Dirac monopole to an SU (2) 
gauge theory defined in five dimensions [20l[2T]. The monopole sits at the origin of a 5D 
space, surrounded by a four sphere S'^. On the S^, this leads to an so(5) homogeneous 
instanton density. In addition to the parameters specifying an ordinary instanton, the 
configuration space of the Yang monopole contains an S"^. This S*^ is the space of possible 
identifications between the gauge theory SU{2) and a space-time SU{2) C 50(5) of the 
isometries of the S*^. A particle charged under this background flux experiences the local 
geometry S*^ x 5*^ geometry. Globally, this S"^ fibers non-trivially over the S*^, so that 
for a single instanton, we obtain a three dimensional complex projective space CPl This 
construction generalizes to SU{n) gauge theory via the principal embedding of the spin 
A^/2 representation of SU(2) in the fundamental of SU{n) where n = A^-|- 1. Via the Ward 
correspondence, the twistor lift of the rank n bundle on S"^, combined with the abelian flux, 
realizes a U{1) x SU{n) gauge bundle Ay on CP^'^. 

Proceeding in the other direction, we can start form the hCS theory with gauge group 
G = U (nNc). Turning on a U{1) x SU (n) gauge field Ay then breaks the gauge group G to 

®See |22] for a discussion of the relation between the 4D QHE and twistor geometry, [23] for 
an embedding of the 4D QHE in string theory and [2lj for a discussion of the connection between 
the QHE effect on and CP'^. For the relation between the lowest Landau level and twistor 
geometry, see also [25] . 
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U{Nc). The low energy limit is a U{Nc) hCS theory, defined on fuzzy twistor space, coupled 
to fundamental matter localized on the twistor lines. From the open string perspective, 
the non-commutativity arises because the open string end points are charged under the 
background flux. The presence of the flux introduces an explicit choice of scale, and the 
open string spectrum is gapped with excitations of energy inversely proportional to this UV 
length scale [S]. The gapped excitations decouple at large N and in the low energy limit, 
the system is forced to lie in the lowest Landau level. 

Our aim in this section will be to study the effective geometry experience by particles 
moving in the lowest Landau level. We begin by reviewing the topology of the Yang 
monopole conflguration. 

3.1 Yang Monopole 

The Yang monopole is a 5 dimensional non-abelian generalization of the Dirac monopole. 
The basic example consists of an abelian flux extended over the full twistor space CP^, 
correlated with a homogeneous SU{2) instanton configuration along the four sphere S^. 
The bundle data is equivalently characterized in terms of the two Hopf fibrations 5^ — CP^ 
and S'^ — )■ S*^ with fibers = U{1) and = SU{2), respectively: 

U{1) -> SU{2) ^ ^7 

; ; . (3.1) 

The basic Yang monopole configuration is given by the original Hopf fibration, which carries 
a single unit of fiux. In this case, the abelian part of the monopole is given by the line bundle 
0(1) ^ CP^ while the non-abelian part is given by a single 50(5) homogeneous instanton 
of SU (2) gauge theory. We refer to the fibrations with fiber 5*^ and as respectively the 
abelian and non-abelian 'parts' of the Yang monopole. 

A concise way to present both contributions is given by working in terms of quaternionic 
coordinates. Given complex coordinates = 'Re{Za) + i\m.{Za) for C^, the quaternionic 
plane can be parameterized as Qi = Z'^ + Z'^j and Q2 = Z^-\-Z^j, where i, j and k denote 
the usual quaternionic generators. Modding out by C* on realizes complex projective 
three-space CP^, while modding out by the left action of the non-zero quaternions on 
realizes the quaternionic projective line HP^ ~ S^. In terms of the Q coordinates, the 
abelian and non-abelian parts of the Yang monopole can now be written as |26] : 

Ay = l {dQlQa - QldQa) . (3.2) 
This represents a U{1) x SU{2) gauge connection via the identification of quaternionic 
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generators with Lie algebra elements of SU{2). More general choices of flux can be accom- 
modated via the principal embedding of SU{2) — SU{n), defined by identifying the spin 
N/2 = [n — l)/2 representation of SU{2) with the fundamental representation of SU{n): 

1 ^ iVl„x„, ^ ^ -2tli, J -2212, k -> -2^13, (3.3) 

where Ifc denote the spin N/2 representation of the SU (2) algebra. Here and in the following, 
the integers and n are related via = n + 1 . The abelian contribution to the Yang 
monopole defines a f/(l) gauge connection 

Auii) = j{dZiZ--ZidZ-), (3.4) 

with curvature equal to times the standard Kahler form on CP'^. The non- abelian 
contribution defines an SU{n) gauge connection over MF^ ~ 5^ which via the Ward cor- 
respondence lifts to an so(5) symmetric rank n holomorphic vector bundle over CP^. 

In the following, the Yang monopole will refer to a gauge connection on twistor space, 
given by the sum of the abelian fiux (13. 4 p and the twistor lift of the SU{n) connection. 
The instanton number of the resulting SU{n) gauge field is maximized by choosing the 
irreducible representation, giving an so(5) homogeneous instanton with instanton number: 

n(n^ — 1) 

"^inst = 7; • (3.5) 
D 

A somewhat more explicit description of the Yang monopole, which also brings out 

its natural relation to the twistor geometry of the four sphere, is as follows |27]. We can 

parametrize the seven sphere S"^ by means of a four-component complex vector (which 
we can view as a spinor of so(5)), subject to the constraint 

Z+Z" = HN, (3.6) 

where for later convenience, we chose the radius of the 5*^ equal to y/hN. For now, however, 
Z and Z^ are just classical coordinates on S'^. The f/(l) action in the first Hopf fibration 
S"^ — )■ CP^ in (13. ip represents the phase rotation 

{Z,Z^) ^ {e'^Z,e-''''Z^). (3.7) 

Similarly, we can parametrize the four sphere with the help of a five-component real 
vector (which we can view as a vector of so(5)), satisfying 

y^VA = f. (3.8) 
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We would like to find a parametrization of in terms of the geometric data of the 
second Hopf fibration S*^ — )■ S*^ in (13.11) . In fact, we have already found this parametrization: 
it is given by the twistor correspondence for 5*^ summarized in section 2.1. Concretely, we 
can package the S'^ coordinates y"^ into a 4 x 4 matrix as in (12. 6p . and impose that Xap 
and Z satisfy the twistor line equation (12. 7p . Indeed, in so(5) spinor notation (and setting 
i = l), the parametrization (I2.12p amounts to setting x = |(1 — T^y^), with the so(5) 
gamma matrices. The twistor line equation (12. 7p then takes the standard form of the Hopf 
equation: Z + y'^V aZ = 0. This equation is explicitly solved, up to an overall phase, via 



V 



Ui 

Ui 
U2 



\ 



J 



HN 



^3 \ni+ in2 



(3.9) 



where the rii are normalized unit vectors on M.^, defining a unit two-sphere n^n* = 1. They 
parametrize the twistor line above the given point y^ on the four sphere. Combined, the 
y^ and rii represent a complete coordinate system on CP^. 

Following [27], we may now characterize the abelian part of the Yang monopole (13. 4p 
as the Berry connection of the spinor Z" defined over CP^. Since the associated holonomy 
is just a phase, the Berry phase Lagrangian for the spinor coordinates is 



-N 



l + ri3 



+ 



v'f.uniy^y'" 

^(^ + ^5) 



(3.10) 



where 77^^ = €ifj,u4^+6i^6u4,—SiuSfj,4 denotes the 't Hooft symbol. Using this, we can reconstruct 
a gauge field flux on both the S*^ direction, and along the S"^ fiber: 



N 



vlu^iy'' 
({i + y^y 



A, = 0, 



A, = N 



(3.11) 



which defines the abelian Yang monopole connection on CP^. 

Both terms in (13.111) have a familiar form. The S"^ part Ai represents the constant 
magnetic flux through a unit two sphere, produced by a Dirac monopole of charge at 
its center. It arises because the two Hopf fibrations in eqn (13. ip are related to each other 
via the basic Hopf fibration that embeds U{1) inside SU{2). Similarly, the 

part A^ looks like the abelianized twistor lift of the basic SU (2) instanton - or rather its 
embedding via the spin N/2 representation inside SU{n) - via the identification 



N ' 



(3.12) 
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Indeed, making this substitution inside of the S'^ part in (13. lip gives a maximal ho- 
mogenous SU{n) instanton on S*^. Again we see that the two components of the monopole, 
the abehan flux (13 ■4p and the non-abelian homogeneous instanton background on S*^, are 
directly linked via the substitution (I3.12p of the 5*^ coordinate with an SU{2) generator 
embedded in SU{n). The complete rank n Yang monopole Ay is the sum of the abelian 
magnetic field (I3.1ip with the twistor lift of the homogeneous SU (n) instanton. 

3.2 Lowest Landau Level 

Let us look a little more closely at the effective geometry experienced charged particles 
in the abelian flux background. At low energy, the particles are forced into their lowest 
Landau level (LLL). This system was analyzed in detail in [9l|27]. The LLL states are 
most succinctly characterized as the Hilbert state obtained by canonical quantization of 
the abelian Berry holonomy Lagrangian (I3.10p . 

We can do this in two different ways. We can first quantize the coordinates, and 
then impose the constraint (13.61) or we can first solve the constraint, as done via the 5"^ x S*^ 
parametrization (13. 9p . Let us first follow the second route. From the second form of the 
Berry action in (I3.10p . we can read off the commutation relations among the Ui and y^, by 
following the usual rules of canonical quantization. One finds coordinates which satisfy the 
following commutator algebra P|[27] 

2 £^ • 2 

[ni,nj] = —eijkUk, [y^.,yu] = ^V^^i, [ni,y^,] = —Vii^uV'' ■ (3.13) 

Rescaling the Ui coordinates that paramterize the two sphere, we see that they turn into 
generators of an su{2) algebra that rotates the S^. Since the rescaled S"^ has radius A^, we 
learn that the su{2) acts via a spin A^/2 representation. The has turned into a non- 
commutative sphere with n = A^ + 1 fuzzy points. Moreover, we see that the su{2) algebra 
does not commute with the position coordinates y'^, but instead act like the generators of 
chiral space-time rotations. The space-time coordinates coordinates also do not commute 
among each other: their commutator is a generator of a chiral su{2) rotation. We will refer 
to the non-abelian commutator algebra (I3.13P later on, when we begin our study of the low 
energy physics of the hCS theory with flux. 

The above discussion also clarifies the mapping (I3.12p between the abelian and non- 
abelian flux: it identifies the spin A^/2 representation of chiral su{2) rotations with su{2) 
gauge generators embedded inside an SU (n) gauge group. The replacement (I3.12p results 
in a homogeneous, so(5) invariant instanton configuration, by virtue of the fact that we 
can combine the chiral su{2) space-time rotations with global su{2) rotations, that acts on 
the Lie algebra labels of the SU{n) gauge field. Indeed one can show that the total space 
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spanned by the Landau level wave functions preserves an overall so(5) symmetry, and that 
moreover, the lowest Landau level transforms as a representation of su{4) [S]. 

We now turn to a more practical description of the LLL states: we first quantize the 
coordinates, and then afterwards impose the constraint fl3.6p . This procedure identifies the 
LLL states with the space of Planck cells, or fuzzy points, on non-commutative CP'^. The 
su{A) symmetry of the LLL level is manifest in this description]^ 

From the first form of the Berry Lagrangian fl3.10l) we immediately read off that the 
four coordinates Z"' act as bosonic oscillators with commutators: 

[Z-,Zl] = h5l (3.14) 

These generate a Fock space obtained by acting with the oscillators on the vacuum 
state |0) annihilated by the Z° oscillators. Each basis state in the Fock space represents 
one Planck cell of the non-commutative space C^, Since this space is non-compact, the 
associated Hilbert space is infinite dimensional. 

Now let us impose the constraint (13.61) . At the quantum level, this will automatically 
lead to a projection from onto the complex projective projective space CP'^. Introduce 
the level operator Hq = Z^^Z"". This operator has an integer spectrum, given by the sum of 
the oscillator levels of the Z-^ number eigenstates. The constraint (13. 6p is now imposed at 
the level of the Hilbert states 

Ho\^) = hN\<il) ; Ho = ZlZ''. (3.15) 

We denote by 'Hcp3(A^) the Hilbert space of states that satisfies this condition. At a 
geometric level, states in 'HQpi{N) represent holomorphic sections of the degree N line 
bundle O^p-i^N). Note that the level constraint Hq = N indeed eliminates one complex 
dimension: it fixes the absolute value of but also implements the U{1) invariance under 
phase rotations — )■ e'"^Z". Eqn (13.151) is the non-commutative way of realizing CP'^ as the 
Kahler quotient C^//f/(l). Observe that the states of 7/cp3(iV) are created by homogeneous 
degree N polynomials in the Z^'s. We can immediately count that 

di„.«,^(iV) = (iV + l)(iV + 2)(iVH-3) ^ 



This formula counts the number of Planck cells that fit inside the compact space CP'^, or 
equivalently, the relative inverse volume of a fuzzy twistor point. We will focus on the 
leading behavior in the limit of large N . In this limit, the size of the Planck cells, i.e. the 
scale of non-commutativity, tends to zero relative to the total size of the projective space 



^The approach to non-commutative geometry we consider has been developed in [28], and we 
refer the interested reader there for additional details (see also [291130]). 
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CP' 



>3|4 



The supersymmetric generalization of the story is straightforward. In addition to the 
bosonic oscillators, introduce four fermionic oscillators satisfying: 



The Hilbert space of points for fuzzy C^'^ is given by the Fock space of states l-L^^m generated 
by z\. The restriction to fuzzy CP^''^ is achieved by introducing the Hamiltonian constraint: 



3.3 Space-Time and Locality 

One of the key features of twistor theory is that space-time physics is a derived notion, 
obtained via the correspondence between complex lines in twistor space and space-time 
points. Since this is a correspondence at the level of the holomorphic geometry, it is 
essentially left intact by the non-commutativity. 

Let us first discuss holomorphic subspaces in CP'^. In commutative geometry, a holo- 
morphic divisor S C CP'^ is specified by the vanishing locus of a degree d holomorphic 
polynomial f{Z°'). In the fuzzy setting, the Hilbert space of points for S is given by states 
of ?^cp3(A^) annihilated by f{Z"'). Let us note that generically this space is non-empty, 
since / : "Hcp^ (^) ~^ ^cp^ — rf) is a linear map to a vector space of smaller dimension. 
Intersections of divisors proceed in a similar fashion. Given two holomorphic polynomials 
/i and /2 of respective degrees di and ^2, the space of states annihilated by both /i and /2 
defines the Hilbert space of points for a fuzzy curve in CP^. Finally, given polynomials /i, 
/2 and /s, we generically obtain a discrete collection of fuzzy points. See [28] for further 
discussion of intersection theory on fuzzy spaces. 

To any point p of commutative twistor space, we can associate a corresponding Hilbert 
space Hp. To see this, note that p can be viewed as the intersection of three linear divisors, 
/^Z" = for z = 1,2,3. The space of states annihilated by all three polynomials defines 
a one-dimensional Hilbert space. This also provides a fiber bundle over the commutative 
CP"^ which we can identify with the abelian part of the Yang monopole. 

Especially important for the twistor correspondence is that for any given point (x"", 6'*") 
on the commutative complexified Minkowski space, there is a corresponding fuzzy CP^(a;, 9), 
with an n = A^-|-l dimensional Hilbert space HxfiiN) spanned by all states \x, 9) that satisfy 
the annihilator equations: 



(3.17) 



(3.18) 



{U^ -ix^'^Tla) \X,9) = 



\x,9) = 



(3.19) 
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These equations implement the commutative twistor hne equations (12.151) as a hnear pro- 
jection on the Hilbert space 'Hx,eiN). This will be a key point of our further discussion: 
although twistor space has become non-commutative, there is still a continuous moduli 
space of twistor lines. Similarly we can define a Hilbert space of dual bra states 'H^g via: 

(X, ^1 {ul - ZXaaTl''^) = ; (x, 9\ {4 - 9,^71^^) = . (3.20) 

The index i of 6ia has been lowered using the infinity twistor X/j. We can thus speak of a 
space of ket states 'Hx,e and bra states "H^^ located at any given space-time point (x, 9). 

Given a specific twistor location (w", tTq, -j/^*) = {x"''^7if,,Ha,9''""iia) on the twistor line at 
{x,9), we can associate an element of Tix^ via the condition: 

(TTa - Atti) A) = 0. (3.21) 

Here A = 7r2/7ri denotes the affine coordinate on the CP^. When combined, the conditions 
(I3.19P and (I3.2ip select a one dimensional Hilbert space within 'Hj.p3|4(A^). 

On the finite radius S*^, there is a closely related way to define position states \y, 9; A), 
where = {y^,y^) denote the S'^ coordinates introduced in eqn (I2.12p . as follows. Begin 
with a state on the CP^I° at the origin where = and y^ = 1. This state |0,0;A) is 
annihilated by u"", and it2 — Xni- Next, we introduce a finite symmetry transformation of 
the supersphere S"''^, which we denote by TZ (y, 9), which maps the origin to the point {y, 9). 
Its bosonic part is an SO{b) rotation. We can thus obtain general position eigenstates by 
applying the symmetry transformation 

|x,^;A) =7^(x,^) |0,0;A). (3.22) 

This defines a unitary parallel transport operation. In other words, each point (y, 9) on the 
sphere S"^'® now comes equipped with an n dimensional linear space l-Lyfi which via (I3.22p 
turns into an su{n) bundle over 5*^'^, or more specifically, the spin N/2 lift of an su{2) 
bundle over S'^. This is the gauge bundle for the non-abelian part of the Yang monopole. 

Since the Hilbert space ^^^^{N) is finite dimensional, it is evident that the states 
associated with different space-time points can not all be independent. Rather, we should 
expect that states at nearby space-time points have a non-zero overlap. The short distance 
scale where locality breaks down, is determined by the overlap between two neighboring 
states \y, 9; A) and {y', 9; \ \. A simple calculation, outlined in [HIH], shows that for y and y' 
close to each other (and close to the south pole where y'^ <C y^) 

{y\ 9, X\y, 9; A) = exp (-^ \y - ) x (0, 9, A|0, 9; A) . (3.23) 
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In the large N limit, the prefactor approaches a regulated delta function, smeared out over 
a small region on S*^'* of linear size 



N 



(3.24) 



As suggested by our notation, this will play the role of the Planck length. 

An equivalent way to see the breakdown of locality is by considering the range of allowed 
angular momenta of functions on S*^. To this end, we introduce a holomorphic position 
operator Xa(3 which acts on a matrix $ as: 

Xo^p^ = Z^^Z\ - Zp^Zl (3.25) 

This definition is motivated by the commutative relation (12. 4p between the space-time 
coordinates Xap and pairs of twistors Za and Wp. The operator Xap indeed selects a 
space-time point, as follows. It is not hard to show [1] that the eigen operators, defined via 

X^p^{x) = x^pX^^{x) (3.26) 

where Xq = I'^^Xap and c-number, are annihilated from the left and right by the 

associated twistor line conditions fl3.19p and fl3.20p . This suggests that we can construct a 
space of functions associated with the space-time as follows. 

Introduce the hermitian conjugate coordinates X^^^^ = Z\a^Zf^]. Starting from the 
vacuum state |0) (0|, we can repeatedly act via X*^. After acting times, we obtain a 
subset of operators on 'H^jpi\4{N), which we denote by Hs-iiN). These matrices form an 
irreducible representation of su{4), given by a Young tableau with two rows of A^ boxes. 
The dimension of this irreducible representation is: 

dim (N) = ^ {N + 3) {N + 2f {N + 1) = h^in + 1) (3.27) 

which is the number of independent functions on 5*^ up to a maximal angular momen- 
tum. We can further decompose this into representations of so(5), which are labelled by 
Young tableaux with two rows of lengths ri and r2 < ri, with dimension (see e.g. |31] ) 
D{ri, = \ + r2 + 2) (ri — r2 + 1) (3 + 2ri) (1 + 2r2) o The functions on a S"^ up to a 
cutoff angular momentum A^ are given by the direct sum of the r-fold symmetric product 
on the vector representation, i.e. representations with r2 = 0: 

N 

Hs^ = ®D{r,0). (3.28) 

r=0 
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Note that LLL states transform in the so(5) irrep with ri = r2 = N/2, with dimension k^- 



21 



In [311132] this space of states was interpreted as the space of functions for a fuzzy S'^. Here, 
we are simply considering a theory on an S*^ where we truncate the spherical harmonics to 
a maximal angular momentum. Since the space of functions scales as N^, the number of 
pixels which can be reconstructed is of order a/ N^. The minimal length scale which can be 
resolved by this angular momentum cutoff is ipi of equation (I3.24p . 

4 Twistor Matrix Model 

After collecting the main geometric ingredients, we are finally ready to introduce the phys- 
ical system of interest. In this section we wish to extract the effective matrix model that 
captures the low energy dynamics of holomorphic Chern-Simons theory in the background 
of the Yang monopole flux. The theory is formulated in terms of finite size matrices, and 
so in particular describes a theory on space-time with a truncated number of degrees of 
freedom. See [33] for earlier discussion on a potential matrix model characterization of 
twistor string theory. 

4.1 Chern-Simons with Flux 

We wish to study the perturbative dynamics of a f/ (nNc) holomorphic Chern-Simons theory 
in the Yang monopole background. As explained, this monopole has two components: an 
abelian flux given by (13. 4p . and a non-abelian flux given by the twistor lift of a homogeneous 
SU{n) instanton with maximal instanton charge kj^f. Both components have an important 
influence on the low energy effective theory. We first discuss the abelian flux. 

The proper way to view the abelian flux is through the lens of the topological B-model. 
In any ordinary Chern-Simons gauge theory, all fields transform in the adjoint of the gauge 
group and the overall U{1) gauge group factor decouples. Hence turning on an abelian 
flux would not have any effect on the low energy dynamics. Instead, we will consider 
the holomorphic Chern-Simons theory as defined via an open topological string theory, or 
equivalently, by taking the zero slope limit of the open twistor string theory introduced 
and studied in [3ll[35]. In the open string theory one can choose to turn on a S-field 
background, which acts like an abelian magnetic field, under which the open string end 
points are oppositely charged. In the zero slope limit, each end point is forced into a lowest 
Landau level orbit, or equivalently, is compelled to occupy a state in '?^^p3|4. 

This gives us a first hint of what the low energy theory should look like: by turning on 
the abelian flux, we have deformed the original U{nNc) holomorphic Chern-Simons theory 
into a non-commutative U (nNc) hCS theory defined on twistor space CP^'^, deformed via 
commutation relations f l3.14p . What does this deformation correspond to in space-time? 
Via the appropriate generalization of the Penrose- Atiyah- Ward correspondence, classical 
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solutions to the non-commutative hCS theory, that is, holomorphic U (nNc) gauge bundles 
with F*^°'^) = 0, correspond to instanton backgrounds on a fuzzy S^. This fuzzy four sphere 
is characterized by the commutator algebra fl3.13p . or equivalently, by the collection of fuzzy 
twistor lines introduced in the previous section. The most symmetric instanton background 
is the non-abelian Yang monopole with rank n and maximal instanton number kN- 

What are the low energy consequences of turning on the non-abelian flux? A first 
obvious consequence is that the unbroken low energy gauge group is reduced from U {nNc) 
to U{Nc). So we should expect the low energy theory to contain a non-commutative hCS 
sector with gauge group U{Nc). But we should also look for other low energy remnants 
of the instantons. From a D-brane perspective, adding an instanton background amounts 
to adding co-dimension four D-branes, one brane per constituent instanton. We will 
refer to the elementary instantons, or co-dimension four D-branes, as defects. Based on our 
experience with intersecting branes, we should thus anticipate the existence of additional 
defect degrees of freedom associated to the non-abelian instantons. An alternative and 
possibly more direct explanation for the existence of defect modes comes from considering 
the ADHM description of instantons. 

One guiding principle for fixing the form of the matrix model is that the low energy 
theory should be compatible with the symmetries preserved by the Yang monopole. Since 
this configuration comes with a choice of scale, specified by the instanton density, we should 
expect conformal symmetry to be broken. However, we do require the bosonic 5*0 (5) 
symmetry to be manifest. Moreover, because of the close link between hCS theory and 
instantons on S"^, we expect that, in an appropriate limit, the configuration space of the 
matrix model should reproduce the ADHM construction of instantons. 

This ADHM correspondence provides important guidance, so let us make the connection 
a bit more precise. Since the Yang monopole breaks conformal symmetry, the low energy 
theory is aware of the radial size i of the S^. There is only one scale because the instanton 
density on S*^ is homogeneous. However, if we give up 5*0 (5) symmetry, we can separate 
the scale defined by the abelian flux from the scale defined by the non-abelian flux. 

In fuzzy twistor terms, this is done as follows. The LLL states admit the action of an 
sm(4) algebra, generated by the hermitian charges J, J, P + K, i{P — K) and iD. The 
hermiticity convention P'^ = K depends on the strength of the abelian flux. The sm(4) 
symmetry is broken to so(5) via the introduction of the infinity twistor (12. 9p . provided we 
make the so(5) invariant choice 7 = 1. As we will see, the appearance infinity twistor is 
linked to the presence of the non-abelian flux. Taking 7 7^ 1 breaks so(6) further down 
to so(4). These infinity twistors are associated with SU{n) instanton backgrounds that 
are homogeneous relative to a four sphere of a different radius, i.e. with non-hermitian 
symmetry generators P + '~fK relative to the inner product set by the abelian flux. Taking 
the limit 7 — )■ corresponds to non-abelian instanton configurations which are centered 
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around the south pole of the S^. In this hmit the abehan flux gets diluted, and we expect 
to make contact with the ADHM construction of SU (n) instantons. This check is performed 
in section [51 

4.2 The Matrix Model 

Let us write the matrix model action. It consists of two terms. The first term is the non- 
commutative holomorphic Chern-Simons action and the second term is the defect action 

>S'mM = S\ics{-^) + Sdefect{Q,Q,A). (4.1) 

^hcs(^) = Tilna^-yD^D^D^^^^ (4.2) 

5defect(Q,Q,^) = Tt(IijQV'QZ'^ (4.3) 

where Qa/B-y is a three index anti-symmetric tensor which is the non-commutative analogue 
of the holomorphic three-form. Here the symbol Tr is the trace over the Hilbert space 
?/j.p3|4(iV) tensored with the u{Nc) color space. We give the precise definition of each ma- 
trix variable in the following subsections. Roughly, each symbol defines a linear operator 
that acts on "Hj^pSK ( A^) ® u{Nc). The = Z°' + A"" are non-commutative versions of 
anti-holomorphic covariant derivatives and is an extension of this derivative to super- 
space. The Q and Q's are the defect modes, and = (Z", ■?/'*) are the non- commutative 
supertwistor coordinates. Finally, X/j is the supersymmetric infinity bi-twistor (12.1 6p . and 
the subscript (•••)^4 indicates the projection onto the top superfield component. 

All matrix variables in the action (l4.ip -f H^ have direct geometric meaning as functions 
and sections of bundles over twistor space. Functions correspond to maps from H^pSK ( A^) 
to 'Hj^p3|4(A^). A convenient representation of such maps is in terms of a power series in the 
oscillators, which we may normal order as M = Mmn\ab{Z'')"^{ip^)°'Z""ilj^. Matrix multipli- 
cation then corresponds to multiplying successive normal ordered power series. Similarly, 
we can describe sections of degree / line bundles by rectangular matrices. These correspond 
to maps from 'H^^3\i{N + p) to 'H^p3\4{N + q), where the net degree is I = q — p. The 
corresponding map has degree A^ + p in the Z oscillators, and A^ + g in the Z"^ oscillators. 
Hence, it can be identified in the commutative geometry with the corresponding section of 
0{l). Note, however, that there is some freedom in how to assign p and g to a degree / 
line bundle. This is because complex conjugation and dualization of a bundle are naturally 
related by Hermitian conjugation of operators. Finally, integration of a product of such 
functions or sections proceeds by tracing over the appropriate Hilbert space |28] . 

In the matrix action, the covariant derivatives D" and V are (0,1) forms on CP^I^ and 
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like the coordinates and (which also act like d operators) and describe non-square 
matrices of degree —1, that change the rank from A^ + ptoA^ + p— 1. The Q's and Q 
variables are sections of a degree 1 line bundle, and thus increase the rank from + p to 
+ p + 1. Modulo such small shifts, all symbols in fl4.ip -f Hl3|l are /cat times /cat matrices. 

4.3 Matrix Chern Simons 

To motivate the form of the holomorphic Chern-Simons action, it is actually helpful to at 
first enlarge the field content, to cover all CP'^'^'s up to level A^. This provides an eight- 
dimensional formulation, as in integral over a ball of square radius ISl inside C^'^, from 
which the hCS contribution arises as a boundary term at level A^. We introduce a fuzzy 
ball B^l^ with a Hilbert space of points: 

-Hb(4|4) = © n„3\,{M) (4.4) 

M=0 

The matrices D" correspond to (0, l)-forms on B^'^. We focus on the bosonic part of this 
differential form content. Each is specified by a power series in the Z'^'s, Z and ip which 
we schematically write as: = D^^^^{Z^y Z"^ {ipy with specified coefficients D"^^^. Let 
us note that this expansion provides a convenient way to deal with matrices which are of 
different sizes. The field content is specified in terms of the direct sum of Km x Km+i 
matrices for M = up to M = A^, which we denote by: 

N 

= .tJ^'^^^-iM+i) (4.5) 

Each term in this direct sum is controlled by the same expansion coefficients Di„^\a- The 
degrees of freedom of the matrix model truncates at finite order set by A^. The transform 
in the adjoint representation of the U {N^) gauge group, and are to be thought of as covariant 
derivatives. To ensure that the gauge fields project to (0, 1) forms on CP^'^, they are subject 
to the constraint [28]: 

ZlD"" = hN. (4.6) 

Expanding around D"" = Z'^ + A^^ this becomes the condition Z]^J^ = 0. This projects 
out the direction orthogonal to CP^'^ inside of C^'"*. See [2S] for further discussion. 

The fuzzy holomorphic Chern-Simons action can now be defined as a trace over the 
8-dimensional ball B^'^ 

^hcs(^) = ^,TT^^^le^,,sD-D^D^D')^^ (4.7) 
where we have introduced g, the gauge coupling of the matrix model. This is the non- 
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commutative version of the continuum action J^^^^ d^^pd'^Zti (J-" A J-') on B"^'^ where J-' = 
OA+AaA. This continuum action is a total derivative, and so integrates to just a boundary 
term, which is the hCS actionH A similar effect occurs in our case: the matrix operator 
inside the trace in (14. 7p naively takes the form of a total commutator: 

e^p^s[D'',D^D^D'] (4.8) 

and thus should have a vanishing trace. However, the D"'s are non-square matrices that 
lower the rank of the Hilbert space by one unit: they map 'Hi^p3\4{N+p) to 'Hf.pa\4{N+p—l). 
So if we cycle around the trace, it lowers the rank of the Hilbert space over which the 
trace is taken on one unit. It is easy to see that, via this mechanism, the trace over all 
states inside the ball B^''^ cancel out, but that one is in fact left with boundary term at 
level A^. This boundary term includes the hCS action as written in (14.21) . 

There is in fact a subtlety in this partial integration argument, that will play an im- 
portant role later on. The eight-dimensional form (14.71) of the holomorphic Chern-Simons 
action is manifestly invariant under the gauge transformations: 

where h = T^Ha is a u{Nc) Lie algebra element and Ha are arbitrary degree zero polyno- 
mials in the and Z's. Hence the result after the partial integration should also be gauge 
invariant. However, as the attentive reader may already have noticed, the action (14. 2 p is 
not automatically gauge invariant, since, even though it is a color singlet, the holomor- 
phic three-form fi^^^ does not commute with general non- commutative gauge parameters 
hA{Z, Z'^). This effect is subleading in and thus disappears in the commutative limit, 
but nonetheless should be taken into account if one insists on writing the hCS action as 
a trace over 'H^f-i\4{N) . There are two ways to deal with this subtlety: (i) stick to the 
covariant action (14. 7p . or (ii) promote VL^fi^ in (14.20 to a dynamical variable that acts like 
a compensator field. We will return to this point in section [71 

4.3.1 Small Phase Space Description 

The twistor matrix model as just given can be viewed as a large phase space description. It 
makes the sm(4|4) symmetry manifest, but at the expense of working with non-square ma- 
trices. There is also a small phase description which sometimes gives a practical alternative. 
Introduce three oscillators for an afiine patch of C^. On this space, we can restrict to 
B^(iV), all states of of level less than or equal to A^. In this description, a homogeneous 
polynomial of degree d in the remaining Z's is replaced by a polynomial of degree d in the 

^^D. Skmner has also considered an eight-dimensional formulation of hCS theory. 
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(^'s. Starting from the bulk action of equation f l4.2p . there is a local direction normal to 
CP'^, which we can identify with Z^. This amounts to the formal replacement of and 
by the identity. Viewing all fields as polynomials in the C's and C^'s, with the maximal 
degree fixed by their homogeneous counterparts, the hCS action simplifies to 

^hcs = ^Tr [E^t^D'D'D"^)^, . (4.10) 
where the trace is over the three ball B'^'^ of states in C'^ with level less or equal to A^. 

4.4 Defect Action 

We now turn to describe the properties of the defect action ( 14. 3p . We first discuss the 
free action, with the coupling to the holomorphic Chern-Simons gauge field turned off. 
So in (14.31) we replace the gauge covariant derivative by the non-commutative twistor 
coordinate . Also, to keep the discussion a bit more transparent, we will for the most 
part restrict our attention to the bosonic sector. The generalization to the supersymmetric 
formulas is straightforward. The free bosonic defect action takes the simple form 

SiQ,Q) = Ti{QDQ), (4.11) 

DQ = I^pZ'^QZ^, (4.12) 

where the trace is over the Hilbert space of points on CP'^ and over the u{Nc) color space. 
This action defines a gaussian matrix model, where the integration variables Q and Q are 
slightly non-square matrices of size ^at+i x /c^v, which in addition may carry a color index, 
transforming respectively in the N^. and N^. of U{Nc). 

The gaussian matrix model (14. lip turns out to have quite remarkable properties, which 
are studied in some detail in the companion paper |T]. Here our main task is to motivate 
why (14. lip is the correct action for the defect modes, that is, for the low energy fluctuations 
around the homogeneous instanton background. So let us list the main characteristics of 
(14. lip , and compare them with our wish list. 

• Global symmetries. Our first requirement is that the defect action should respect the 
5*0(5) symmetry of the homogeneous instanton background. The action (14. lip clearly does 
- but it is instructive to see how it works out. Let us first look at the global symmetries of 
the kinetic operator D defined in (I4.12p . The operator D is in fact invariant under a full 
set of gl{4:, C) transformations. Consider a gl{4:, C) generator A^^/j = Zj^Z/s, and define its 
action on Q via M-apQ = -^apQ — Q-M/Sa- Note that the indices a and f3 switch locations 
between the two terms. A short computation shows that [D, A^°^] = 0, and so the kinetic 
operator is invariant under the group generated by all A^°« operators. It would not be 
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correct, however, to conclude that the defect action preserves (7/(4, C): the action (14. lip 
involves a trace, which is defined with respect to a choice of inner product. So the most we 
could hope for is the symmetry algebra m(4). Additionally, the action contains the infinity 
bi-twistor lap, which transforms as a vector of so(6) ^ su{4). This breaks the symmetry 
group to so(5). Concretely, hermitian generators that leave the action invariant are 

M{v) =v''pZlZ^, ^;"/3 = y/3«, (4.13) 

where the raising of the second index is accomplished with the J'^'^ bi-twistor. This leaves 
10 hermitian charges, which are the symmetry generators of so(5). In a similar way, we 
find that the supersymmetric action (14. 3 p preserves the symmetries of the supersphere S^^^. 

• Ultra-locality. A second desired property is that the defect action should refiect the 
quantum Hall intuition, that its excitations are bound to Landau orbits localized along the 
5*^. Alternatively, using the twistor string theory terminology of [7], we wish to see that 
the defect action shares the properties of an effective action for (the ground states of) open 
strings that stretch between the instantonic D-branes and the space-filling D-brane. Both 
arguments indicate that the defect action should be ultra-local in space-time. Eqn (14. lip 
satisfies this beautifully [T]. Consider the position operator X^js introduced in eqn (I3.25p . 
An easy computation shows that these operators commute with the kinetic operator: 

[D,Xap] = 0. (4.14) 

The same is true for the hermitian conjugate operators X^^. Hence the kinetic operator 
acts within the eigen space (I3.26P of Xap. It is the precise sense in which the defect action 
is ultra-local along the S*^. So it has the right to be interpreted as the action of modes that 
are tied to local defects, that wrap the twistor lines. It is easy to verify that this result 
extends to the supersymmetric case. 

• CP^ Propagation. This ultra-local property indicates that the defect modes propagate 
along the CP^ fiber directions. Continuing to take guidance from twistor string theory, we 
would like to see that, when restricted to a particular fiber, the action reduces to that of 
a 2-dimensional chiral free field. Via 5*0(5) symmetry, it is sufficient to verify this for the 
twistor line at the origin. So let $ be a field localized at the origin, i.e. it is a linear map 
on "Hj-psK that is made up only from tTq and vr^ oscillators. On this subspace, the kinetic 
operator D reduces to 

D$|„i=A,<l>7rb (4.15) 

Via the commutation relation [tTq, tt^] = Sab, we immediately see the right-hand side indeed 
acts on $ as a standard Dolbeault operator d = vr^Tp^- on CP . Based on this alone, we 
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expect that the propagator of the D operator looks hke the propagator of a two-dimensional 
chiral free field. The propagator is studied in some detail the companion paper pQ, where 
this expectation is confirmed. For additional discussion of the Dolbeault operator on fuzzy 
CP\ see e.g. [36H3H]. 

• ADHM correspondence. The defect modes Q and Q have a familiar analogue the 
ADHM construction of instantons. From the comparison with ADHM we learn that they 
naturally transform in the fundamental and anti-fundamental of u{Nc). In the full system 
described by the action fl4.ip -( 10|) . they provide a source term to the equation of motion 
of the holomorphic Chern-Simons gauge field. Looking at f l4.3p and the small phase space 
form f l4.10p of the hCS action, we can already see that the (top component of) the D°' 
equation of motion takes the schematic form of the ADHM equation. We will make this 
match more precise in the next section. 

• MHV correspondence. The gauged defect action is invariant under gauge rotations 

Q e'^Q, Q ge"^^ ^ e^^P^e"^'^. (4.16) 

We read off that the bulk gauge field couples to the current: 

= Tf^Q'ZpQ^ (4.17) 

where is a U{Nc) generator in the fundamental representation. This current is conserved 
on-shell: Z^J^ = 0. Hence, after taking the large limit, the matrix model has all elements 
in place to provide a candidate dual description of 4D Af = 4 scattering amplitudes. In 
particular, following [39] and [7], we are led to identify the correlation functions of these 
current operators with MHV gluon amplitudes. As shown in P, the two sides indeed match. 
We will study this relationship in more detail in section 6. 

We have verified that the defect action fl4.3p passes several non-trivial checks, which 
support its candidacy as the correct action for low energy fluctuations around the instanton 
background. In the next sections, we will study the physical properties of the matrix model 
in more detail. In particular, we will present evidence that, in a suitable large A^ limit, its 
correlation functions approach those of A/" = 4 gauge theory, coupled to gravity. Various 
technical details will be delegated to p. 

5 ADHM Limit 

Let us briefly review the realization of ADHM in terms of bound states of branes. The brane 
construction of ADHM is discussed in [10Hl2] (see [13] for a concise review). To frame our 
discussion, consider type IIB string theory on M^'^ with n D7-branes filling R'^'^ x M^, and k 
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probe D3-branes filling M^'^. We focus on the directions, relative to which the D3-branes 
look like instantons. 

The theory on the probe D3-branes is given hj a. U{k) Af = 2 supersymmetric field 
theory. In A/" = 1 language, the motion of the D3-branes is parametrized by three chiral 
superfields, (pi, 02 point parallel and ip points normal to the seven-brane. All three fields 
are in the adjoint of U{k). Additionally, there are bifundamentals q and q which define 
n X k and k x n matrices. These correspond to 3 — 7 strings. The superpotential reads: 



In this system, the trace over U{n) is associated with the flavor symmetry of the D7-branes. 
There is also a D-term potential. To study the classical vacua, it is enough to consider 
solutions to the F-term equations, modulo the complexified gauge group GL{k,C): 



These are the ADHM equations. The moduli space consists of several branches: the 
Coulomb branch, the Higgs branch, and mixed branches. The Coulomb branch, q = q = 0, 
describes the motion of the D3-branes away from the D7s. On the Higgs branch, (f = 0, 
and q and q are non-zero. In this case, the D3-branes have dissolved as instantons inside 
the D7-branes. On the mixed branches, some combinations of both are non-zero. 

The Higgs branch coincides with the A;-instanton moduli space of the GL{n, C) gauge 
theory on W^. This case is characterized by [01,02] = qq, and the space of solutions 
modulo GL{k, C) has real dimension 4kn, which is the dimension of the instanton moduli 
space@ More generally, there can be isolated solutions with all fields switched on. From 
the perspective of the D3-branes, such mixed branches describe configurations where the 
D3-brane has partially puffed up in the directions normal to the seven-brane. 

Let us now see how this matches up with the behavior of the matrix model. To do this, 
we will consider a slightly broader notion of the system, where all modes are complexified. 
We can see roughly how the match works. There are three independent adjoint fields 0i, 02 
and (f, which are naturally identified as the covariant derivatives Z)". In a local patch 
of the 5*^, we can pick a preferred origin and a corresponding M^. Equipping this with a 
complex structure we have a C^, and two corresponding vector fields Di and D2. The third 
component Dq points along the direction of the twistor line at the origin. In addition, the 
defect modes Q and Q have their analogues with the D3-D7 bi-fundamentals q and q. Note, 

-"^^The counting is as follows. We have 4fc^ + 4/cn real degrees of freedom from 0i © 02 and q® q. 
The F-term constraints impose 2k'^ conditions, and modding out by GL{k, C) removes a further 
redundancy of 2k'^. The resulting moduli space has dimension Akn. On S^, the instanton moduli 
space is 4/cn — n^, but if one allows gauge transformations at infinity, the dimension is again 4/cn. 



Wprobe = Trc/(fc) {(f [01, 02]) - Tluin) (q'fq) 



(5.1) 



[01, 02] = qq, l^P, 02] = [v, 0i] = 0, q(p = (pq = 0. 



(5.2) 
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however, that the size of the matrices do not yet match: Q and Q are k^+i x matrices, 
whereas q and q are k x n matrices. 

We now provide a more precise match with ADHM. First let us specify the hmit in 
which we should expect that match to become exact. For this we need to be able to turn 
off the non-commutativity. This is done as follows. Recall the form of the infinity bitwistor: 

^iJ = l\ e'^' (5.3) 
V r],J 

Let us consider the limit 7 — )■ with i held fixed. As explained, the 7 parameter describes 
the ratio of two radii, namely (i) the radius of the S"^ relative to which the abelian flux 
of the Yang monopole is homogeneous, and (ii) the radius of the S*^ relative to which the 
SU{n) instanton part is homogeneous. So sending 7 — )■ amounts to localizing the abelian 
flux in a small region near the north pole relative to the non-abelian flux, or equivalently, 
localizing the non-abelian flux in a small region near the south pole relative to the abelian 
flux. From either perspective, it is clear that in this limit, the effect of the abelian flux 
on the non-abelian configuration becomes negligible. Let us take the perspective where we 
zoom in on the south pole region of the S*^. In this limit, the four sphere starts to look like 
fiat space time, and moreover, only the non-abelian part of the configuration is retained. 
Varying the matrix action (14. 3 p with respect to Q, we obtain the equation of motion 

VaQ = TijV'QZ' = (5.4) 

When Xjj is invertible, this equation generically possesses no solutions. Indeed, via the 
rank-nullity theorem, a map : /'^at+i x Kj^j — > x J'^at+i will in general have a trivial 
kernel. In the specific limit 7 — > 0, however, Xjj is no longer invertible and the matrix 
system will have a moduli space of solutions. Looking at eqns (15. 3p and (15. 4 p we see that 
derivatives in the uj direction cost much less energy than derivatives in the other directions 
of the supertwistor space. The effective geometry experienced by the modes is C^'^ x CP^'° 
where the CP^'" factor is composed from the oj directions. (In the presence of a background 
gauge field, the overall direction of this constraint will change, but this can be absorbed 
into the definition of the oscillators.) Without loss of generality, the constraint is then: 

TTaQ = ^'Q = (5.5) 

These conditions project Q down to a (ra + 1) x Kj^ supermatrix, which is purely bosonic 
in the {n + l)-component part. In other words, the left index on Q are states of a bosonic 
CP^'° at level n. The kinetic term for these modes then fixes the form of Q, so that it is 
given by a -R'at+i x n supermatrix. 
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In this limit, it is also appropriate to take a specific basis for the bulk gauge fields. 
There are the two derivatives in directions transverse to the CP^'*^, given by and D^. 
In addition, there is the derivative along the CP^'°, given by D^. To analyze the remaining 
equations of motion, it is helpful to pass to a description in terms of an affine patch. This 
leaves us with three matrices D^, D^, which admit a T/^-expansion in terms of x kjsf 
bosonic matrices, and Q and Q which respectively admit ip- and ip"^ expansions in terms of 
bosonic matrices of respective sizes {n + 1) x kj^j and kj\f x (n + 1). The resulting action is: 

S = Tr,^(D° [D\D'])^, + Tt^{QD'^Q) (5.6) 

The ADHM equations now naturally appear by varying with respect to the tp'^ compo- 
nent of = + . First, observe that ip'^W is sandwiched between the supermatrix 
parts of Q and Q. Since ip"^ corresponds to an annihilator operator when acting to the right 
on Q, the only non-zero coupling involves the ip'^'^ component of Q, which we denote by q. 
On the other side, the only term that survives is q, the purely bosonic component of Q. 
The terms of the bosonic action involving 6° are then: 

S = Tr„3 (feO [d\ d^] ) + Tr„i|„ (5.7) 

where d^ is the bosonic component of D^. Setting to zero all other components of the 
supermatrices, we see that the resulting equations of motion are identical to those of line 
( 15. 2p . Working modulo GL{ki^, C) transformations, we obtain a moduli space of instantons 
for su{n + l) gauge theory at instanton number /c^r. Note that this is the maximal instanton 
number for which it is possible to define an so(5) homogeneous instanton configuration: 

(iV + l)(iV-f2)(iV + 3) {n + l){{n + lf-l) 
kn - = ^ • lo-«J 

This is again in accord with the fact that we are considering instanton configurations which 
are deformations of the non-abelian part of the Yang monopole. 



6 Continuum Limit 

Our discussion so far has focussed on the finite regime of the matrix model. The finite 
matrices furnish a basis of functions on 5"^ but only up to some cutoff angular momentum. 
To arrive at a continuum theory, it is necessary to take a large limit. 

In addition to TV, the matrix model comes with a length scale, given by the S'^ radius 
and a parameter 7 which controls the breaking of the so(5) symmetry. Depending on 
the scaling of these parameters, we arrive at different limits for 4D physics. The most 
basic large limit is taken while keeping 7 = 1, which preserves the so(5) symmetry. 
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Alternatively, we can take a combined large N and fiat space limit, which corresponds to 
performing a Wigner-Inonii contraction of the original so(5) algebra. This involves rescaling 
the momentum and conformal boost generators, or equivalently, rescaling the tt''" oscillator 
relative to the oscillator. 

We can consider two types of fiat space continuum limits. In the first scenario, we 
first send to infinity, and then take the fiat space limit i — oo. As we will argue, 
this continuum limit produces pure A/" = 4 SYM theory. The correspondence arises in a 
similar way as in twistor string theory, and the dictionary looks almost identical. But there 
are some important contrasts, since our treatment of the defect system is quite different 
from that in [7]. Most significantly, we will present clear evidence that in our set-up, the 
continuum A/" = 4 gauge theory arises without any coupling to conformal gravity. 

A second possibility is to take a double scaling limit, in which both A^ and i are sent to 
infinity, but such that the ratio = £'^/N is held fixed: 

N ^ oo,i ^ OO, with = — fixed . (6.1) 

We will argue that in this limit, the system does contain gravity. Moreover, we will present 
evidence that the emerging gravitational theory is described by the usual Einstein action, 
with a finite Newton constant Gat of order i^i. The strongest evidence in support of this 
inference is that the matrix model is able to reproduce the MHV graviton amplitudes. The 
calculation that leads to this result is described in the companion paper [T]. In the next 
section we will try to give a more conceptual explanation of how and why the gravitational 
degrees of freedom can arise from a matrix model. 



6.1 Emergent Gauge Theory 

Let us first begin with a discussion of the purely gauge theory sector. In the matrix model, 
this is arranged by taking ipi — t- 0. Here we will study this limit directly at the level of the 
action. 

The basic dictionary is relatively simple and standard. Starting from the matrix model 

action fl4.1l) - fl4.3p . the first step is to rescale all twistor coordinates by a factor of so 

that the commutation relations take the form [Z", Z^] = S'j^/N, etc. Hence, upon taking the 

large A^ limit, the twistor coordinates become ordinary commutative variables. Accordingly, 

the trace of the Hilbert space 'H^^iii of points, reduces to an integral over the commutative 

supertwistor space CP^'^. The covariant derivatives D", instead must be defined with an 

d 

extra prefactor of A^, so that reduces to a anti-holomorphic derivative We will 
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thus write 



where A°' denotes the gauge field. Applying this dictionary to the matrix version of the 
holomorphic Chern-Simons action gives the continuum hCS action, in the form explained 
at the end of subsection 14.31 In particular, the stationary points of the hCS action are 
holomorphic fiat connections which satisfy fl2.19[ 

The appearance of a continuum field theory in the large N limit, while standard and 
reasonable, in fact needs a bit more justification. For any finite N, regardless how large, 
the space of functions on a non- commutative space includes arbitrarily non-local maps. 
Indeed, the matrix model involves integration over matrices D", Q and Q, which include 
maps from one point to any other point on the fuzzy CP^'^. In our context, these non- 
local maps are exponentially suppressed, due to the specific form of the matrix action. 
For example, as seen from f l2.19p . the saddle points of the matrix hCS action are maps 
that are locally holomorphic in the Z"'s0 Since the Z°"s commute among each other, 
holomorphic functions multiply locally. In other words, all non-locality is associated with 
the dependence of the fields. As we will see more explicitly in the next section, the size 
of this non- locality is set by ipi. So in the ipi — )■ limit, the theory becomes local. 

The correspondence with local 4D physics proceeds by projecting the CP'^'^ gauge field 
A'^ down onto S^. This is done via the Penrose correspondence, which repackages the 
components of A"' into those of A/" = 4 SYM theory. Adjusted to our context, the Penrose 
transform works as foUowsH 

Suppose A is an on-shell mode. So on a suitable subset of the Hilbert space, it commutes 
with the coordinates. In particular, it preserves the twistor line equations f l2.15p that 
select the -|- 1 dimensional Hilbert space Hx associated with a given space-time location 
X. (Here and below x is short-hand for a point {x,9) on the supersphere S"^'®.) We can 

■^Functions or sections of degree p line bundles on CP^I^ can typically not be globally holomor- 
phic, since it should be expandable as a degree zero (or degree p) polynomial in terms of and 
zt- A function f{Z) locally holomorphic is if its commutator with the holomorphic coordinates 
is an operator that annihilates a subspace 'H{S) of the Hilbert space T-L^p3{N). This subspace 
'H{S) represents the local region within which we can consider f{Z) as a holomorphic function of 
the Z coordinates. For a more precise definition of locally holomorphic functions, see [T] 

^^In commutative twistor theory, the basic correspondence is between space-time fields of helicity 
h, and certain cohomology elements on PT', that is, CP^ with the CP^ at infinity removed: 

{helicity h space-time fields} = (PT', 0{2h - 2)) (6.3) 

where this is defined with respect to a sheaf cohomology. Consider a meromorphic (0, 1) form 
defined over commutative twistor space, which we denote by ^{u;,Tr). Imposing the constraint 
_ ix'^'^TTa, we obtain a space-time dependent (0,1) form ^'(ix7r,7r). Given this, we obtain a 
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therefore decompose ^ as a sum of operators 




(6.6) 



where ^"(x) denotes a hnear map from Hx to Hx- As motivated in more detail in [T], an 
on-shell mode corresponds to taking ^"(x) of the special form 



where, in order for A to represent a holomorphic function of the twistor coordinates (a;", iTa), 
A must be a function of the form A{ix°''^Xa, Xa)- Equation (16. 7p is the quantum version of 
the Penrose transform. The classical version is obtained by sandwiching both sides between 
the two position eigenstates (x, 0| and |x, 0), and using that (x, 0|x, A) = (x, A|x, 0) = 1. 

Applying this reduction procedure to the pure hCS theory only gives the self-dual gauge 
theory on S*^. The role of the defect modes is to supplement the rest of the A/" = 4 SYM 
theory. At the level of perturbation theory, the match proceeds via the CSW rules for 
constructing scattering amplitudes [H]. To see how this emerges from the large matrix 
model, let us look more closely at the defect contribution. 

As we have seen, the defect action is ultra-local in space-time. The fields Q and Q do 
not propagate along the space-time directions, and hence do not correspond ordinary space- 
time fields. The defect modes appear quadratically in (14. 3p . So it is natural to integrate 
them out, and collect their contribution to the effective action for the space-time gauge 
field A in the form of a functional (or rather, for us, just a regular) determinant. 

Let us write the defect action (14. 3 p in the short hand notation 



with D the free kinetic operator defined in eqn (I4.12p . The trace is over the Hilbert space 
7i^p3\4,{N) of point on twistor space. The (lowest superfield components of the) defect modes 
Q and Q are both k^+i x matrices. In the following, the space of all such matrices will 
be denoted by M^. Now consider the operator D oD_4 = 1 + D A. This defines a map 

helicity +n or —n free field on ordinary space-time via the Penrose transform: 




(6.7) 



'S'defect — 



Da 



D + A 



(6.8) 




(6.4) 



(6.5) 



where the integration is over a contour 7 defined over the CP 



in the vr coordinates. 
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from Mk to M^. So we can consider its determinant. Integrating out the defect modes leads 
to an ^-dependent effective action 



= Tr log(l + D k), (6.9) 

where we used the standard identity log det = Tr log. We will now rewrite this in the form 
of a generating function of MHV gluon amplitudes. 

Let us specialize to the case that A describes a sum of on-shell modes, of the form (16. 6p - 
(16. 7p . The defect mode kinetic operator D and the on-shell field A then both commute 
with the position operators Xq,^ 

\D,X^p\=Q ■ [A,X^p]=Q. (6.10) 

The operator inside the trace in (16. 9p also commutes with X^^, and thus defines an ultra- 
local operator on space-time. Given this result, we can naturally decompose the trace over 
Mfc into two factors: a trace over the space-time directions, and a trace Tr^ over the space 
of maps that act on the N +1 dimensional Hilbert space l-Lx associated with the space-time 
point X. In the large N limit, this decomposition amounts to making the replacement 

Ttm^ ^ I d^^^x Tr^. (6.11) 



Here we have made use of the fact that the overlap between states located at different 
space-time points takes the form of a sharply peaked gaussian (I3.23p . which in the large 
limit turns into a delta-function. In the continuum theory, (I6.1ip simply amounts to 
writing an integral over CP^'^ as an integral over the 5*^'^ base, times an integral over the 
CP^ fiber. The D operator in (I6.13P takes the form (I4.15p . which in the large N limit 
becomes the standard Dolbeault operator d on CP^ In the companion paper P, we work 
out the precise form of the Green's function D and show that at large TV it coincides with 
the continuum version. We can view D as the kinetic operator of a free chiral field theory 
in two dimensions. To reflect this, we shall make the replacement D = d when acting on a 
given twistor line. 

We thus conclude that (16. 9p in fact describes a local space-time action 

Sos{A) = j(f\\C,^{A{x)) (6.12) 
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The effective Lagrangian takes the form 



C,s{A{x)) = Tr,log(l+9 



(6.13) 



logdet^.(l + c> A{x 



)) 



(6.14) 



where Tr^^ and det^; respectively denote the trace and determinant over the space of hnear 
maps on Tix, or in more geometric terms, the space of functions on the commutative hmit 
of a fuzzy CP^'° twistor hne located at x. Equation fl6.13p is then the generating function 
IPI correlation functions of bi-linear currents of the chiral free fields, of the form f l4.17p . 
In [T] we show that the large N limit of these current correlation functions reproduces the 
Parke- Taylor formula. 

At this point we have made contact with the works P^HB]. where it is shown that the 
continuum version of the effective action fl6.12p -( l6.13p is the generating function of gluon 
MHV amplitudes. Importantly, as pointed out in these references, the effective action 
(I6.13P evades the troubhng appearance of conformal gravity amplitudes, that plague the 
twistor string proposal of [71117]. The main difference between the two proposals is that 
in the twistor string theory, the effective action that gave the MHV amplitudes arose from 
integrating out the open string modes attached to defect Dl-instantons, which then were 
integrated over their moduli space of positions. This latter procedure gives an effective 
lagrangian of the form det(l + d A), rather than logdet(l + d A). The det expression 
of twistor string theory has several drawbacks. First, it is not by itself gauge invariant, and 
to be well-defined, requires additional couplings to a closed string sector. By contrast, the 
log det term is gauge invariant. In the continuum theory, this is due to a subtle interplay 
with A/^ = 4 superspace0 For us, the gauge invariance of fl6.13p is manifest, since it was 
derived by integrating out a regulated and manifestly gauge invariant action fl6.8p . 

Let us note that when there is reduced supersymmetry, the cancellation of chiral anoma- 
lies is a bit more subtle. In this case, the 2D chiral determinant in (16.131) would produce an 
anomaly outflow, that somehow needs to get cancelled by some anomaly sink somewhere 
else. The appearance of the anomaly is related to subtleties with the analogue of equation 
(16. lip . With less than J\f = A supersymmetry, the CP^ Hilbert space does not neatly fac- 
torize as assumed in (16. lip . The obstruction is the Berry phase associated with passing 
from one coherent state \x) to another \x'). The associated Berry connection is the Yang 
monopole background. The Yang monopole background transforms under gauge transfor- 
mations, and in this way, is capable of transporting the anomaly from one twistor line to 
another. So we see that the matrix model regulates the log det action proposed in [l5lll6] , 

^^Indeed, without the superspace integrations, the MHV action of |45p46j would not have been 
gauge invariant, and moreover, would have contained additional divergences due to ultra-local 
interactions and accompanying factors of S{0). 
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and should make sense with reduced supersymmetry. 

Returning to the A/" = 4 symmetric case, one could perhaps wonder how it could be 
possible that in the Ipi — )■ limit, one can recover a superconformally invariant continuum 
theory. At any finite only the 5*0(5) symmetry is unitarily realized in the matrix model, 
so where could the other generators come from? The mechanism is expected to be the same 
as for any non-conformal theory with a (non-trivial) IR fixed point: conformal symmetry 
is broken as long as the UV scale is present, and arises as a linearly and unitarily realized 
symmetry only in the strict IR limit. 

Finally let us comment on non-MHV amplitudes. The basic story here is the same as 
in the continuum theory. The total matrix action combines the holomorphic Chern-Simons 
action with the effective action fl6.12p - fl6.13p . As in [HlllH], general types of amplitudes are 
obtained by gluing together MHV amplitudes via propagators of the hCS theory. Here we 
do not have much to add to the continuum discussion, except for the folUowing comment: 
instead of integrating out the defect modes, we can alternatively integrate out the bulk 
modes. Specifically, one can choose an axial gaugj^. for which the interaction terms of the 
hCS theory drop out. So we can just perform the gaussian integral over the gauge field. 
The effective action for the defect field Q and Q will then contain the kinetic term, and 
a non-trivial quartic interaction term, which can be included iteratively. It is tempting to 
conjecture that the BCFW recursion relations |19] can be interpreted as the loop equation of 
this interacting matrix model. We leave the exploration of non-MHV and loop amplitudes 
for future study. 

7 Emergent Gravity 

In the previous sections we have provided evidence that the continuum limit of the matrix 
model specifies a U{Nc) gauge theory on an S^. We have also seen signs that hint at the 
presence of gravity. This is refiected in the appearance of a minimal length scale ipi, and the 
fact that gauge transformations on the non-commutative space look like diffeomorphisms. 

In this section we assemble these hints into concrete evidence that the matrix model 
contains a gravitational subsector. In the companion paper [T] we show that MHV graviton 
scattering amplitudes can be represented as correlation function of currents of the gaussian 
matrix model, that describes the defect modes in a fixed A background. In this section, we 
will try to embed this result inside a more global geometric perspective. The appearance 
of a gravitational subsector is in some sense anticipated by earlier work on twistor string 
theory. However, as opposed to there, the matrix model comes with an intrinsic length 
scale, and therefore does not seem to describe conformal gravity. We comment on this 

Another natural and potentially useful gauge condition on A is Z'^Aa = 0. Note that this 
condition involves the infinity twistor. 
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point further in subsection 17.41 In this section, we will put Nc = 1. 



7.1 Currents and Compensators 

In this subsection we show that the matrix model contains gravitational current operators. 
To some extent, the presence of this sector is necessary in order for the non- commutative 
Chern-Simons action to even be well-defined. As explained earlier, starting from the hCS 
action (14. 7p . written as an integral over the four-ball B^'^, one could in principle perform 
a partial integration and obtain an action of the form (14. 2p . However, the naive action 
(14.21) is not invariant under the gauge transformations D°' — e'^'^D°'e~^'^, since the /I's 
do not commute with Qa/s-y- To rectify this problem, we can introduce a compensator 
field = ^e'^^^^QiBys which transforms in the adjoint under non-commutative U{1) gauge 
transformations. The correct action is then: 

5fl.cs = \e^fs^sTT {n-D^D^D') . (7.1) 

g 

One can view the fi" as a residual boundary or 'singlet' component of the U (Nc) gauge field 
on the eight dimensional ball B^'^, which is left over after performing the partial integration. 

This gives a clear indication that the adjoint f/(l) factor of the fuzzy system is on a 
different footing from the other gauge symmetries. In commutative terms, maintaining 
gauge invariance via the compensator can be interpreted as allowing the holomorphic form 
fl = eai3-y5Z°'dZ^dZ'^dZ^ to be deformed, akin to a Kodaira-Spencer theory of gravity [50] . 
In so doing, the geometric data of the twistor space become dynamical. This already 
suggests the appearance of a gravitational subsector. 

It is reasonable to assume that the compensator should also enter in the defect action. 
Indeed, there is a natural place for it. Viewing the defect action as descending from a 
holomorphic Chern-Simons action in the presence of a flux defect (as is common in the 
study of intersecting brane configurations), the compensator appears in the action as: 

Saefect = Tr (j„^gD°gfi^) . (7.2) 

The Q and Q modes have become bifundamentals of u{Nc) x u{l). The system is invariant 
under the enlarged group of gauge transformations: 

Q^e'''Qe~'^, Q^e'^Qe-''', ^ e*'^D"e-*\ fi" ^ e*^^fi"e-*^^. (7.3) 

where his a. Lie algebra element of u{Nc) and g is a phase, and both are local functions of Z 
and Z"^. In addition to the usual color gauge rotations, we see that the gauge transformations 
can also shift the locations of fuzzy points. In other words, the full symmetry is u{Nc) x 
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n(l) X gl{kN) X gl{kN) where the two gl{k) factors act by left and right multiphcation 
on the defect fields. In keeping with the interpretation of the fi" mode as a compensator 
field, it is natural to focus on the diagonal gl{kN)diag C gl{k]^) x gl{k]\f). Indeed, if we 
view the matrix Q as a field located at some holomorphic point x, then the simultaneous 
left and right action by the diagonal gl{kN)diag corresponds to a diffeomorphism on the 4D 
space-time which is capable of shifting x to another location. 

Following the Noether procedure, we obtain the gl{kj^) x gl{kp^) 'currents' (here we 
switched notation, and explicitly write the glik^) matrix incdices) 

MT) = r^iQ''~'ZpQ'^\ Mf) = tiQ'^-^ZpQ'^K (7.4) 

where T^i is a generator of gl{k]\f). The diagonal currents are the linear combinations 

^^{T) = MT) - MT). (7.5) 

In principle, one can consider charges associated with arbitrary glikjsf) transformations, in- 
volving an arbitrary number of oscillators. However, we can anticipate that the equations 
of motion of the corresponding bulk gauge field disfavors such anti-holomorphic dependence. 
So it is reasonable to restrict our attention to currents of the form 07.41) which contain only 
one single Z"^ oscillator. We will write out these currents in the next subsection. 

Associated with the diagonal gl^kj^j) is a gauge field 01^ , which is a function of Z and 
Z"!". It couples to the defect system via the vertex operator: 

r(2t) = Tr (21^ J^j) - TiiOl^Jp) (7.6) 

Based on the index structure of the diagonal glik^) charges 07.41) . we see that the 
symmetries define (0, l)-forms valued in the holomorphic cotangent and tangent bundle 
to twistor space. This is the same mode content expected for a conformal graviton |47j . 
but with a few important differences. First, we see that they couple to a defect action 
which involves the infinity bitwistor. This means that although these currents are obtained 
from conformal gravitons, they are subject to the constraint that they leave the defect 
action, and thus the infinity twistor invariant. This imposes conditions which truncate the 
physical mode content. The combined presence of a UV cut-off and of the infinity twistor 
both break conformal invariance, which suggests that the emerging gravitational sector will 
not be conformal. 

We can extend our discussion to the supersymmetric case with defect action: 

S-defect = Tr {XijQV'Qn'^ . (7.7) 
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The bosonic currents (i.e. with / and J bosonic) provide us with the gauged translations 
and u{Nc) rotations. The currents with I and J both fermionic give rise to su{A) gauge 
fields. Further, there is a fermionic current as well which couples to the gravitinos. All this 
is in accord with the gauge symmetries of an A/" = 4 supergravity theory!^ 



7.2 MHV Graviton Scattering 

To confirm the proposed geometric interpretation of the gl{k]sf) symmetry, in the companion 
paper [T] we show that correlation functions of the currents fl7.4p reproduce MHV graviton 
scattering amplitudes. Here we just briefiy comment on how this structure emerges from 
the matrix model, and on how it relates to existing results obtained from ordinary twistor 
space. 

MHV amplitudes in gravity involve a single ingoing negative helicity graviton and an ar- 
bitrary number of ingoing positive helicity gravitons. Collectively, the plus helicity gravitons 
represent a selfdual background geometry, on which the minus helicity graviton propagates. 
The MHV amplitudes arise by expanding out the background field in terms of linearized 
perturbations around fiat space-time. The twistor realization of this calculation builds on 
Penrose's non-linear graviton [18]. In this construction, (anti-)selfdual space time back- 
grounds are represented as complex structure deformations of twistor space. Infinitesimal 
deformations correspond to vector fields, that via their Lie derivative, bend the location of 
the holomorphic twistor lines. As shown in [51], this geometric description allows a com- 
putation of MHV graviton amplitudes, that reproduces the well-known BGK formula |52j . 

The basic setup for scattering theory in the matrix model is to zoom in on a small and 
locally fiat patch near the south pole of the S*^. This allows us to define a fiat space limit 
via the Wigner-Inonii contraction of the so(5) algebra: 

V^P + r'^K (7.8) 

so that locally, translations are given by just the P. As shown in [T], the asymptotic data of 
the scattering theory are specified by vertex operators T such that the adjoint U{1) gauge 
field takes the special form 21^ = Z^T, where T{Z, Z^) is some non-commutative function. 
The gl{kN) generators of interest for MHV graviton scattering are of the form which locally 
have a single oscillator: 

T„(p) = (i;-P)^(p) (7.9) 
where v'^°' is a polarization tensor for the graviton, v-V = v"'°'Paa + VaaK"'"' is a complexified 

^^The bulk gauge field also can be viewed as a vector on superspace. We are making 
the simplifying assumption that these additional components of the would be gauge field in the 
fermionic directions have been set to zero. It would be interesting to study this more general 
situation, though it is unclear to us that it is a necessary element of the construction. 
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so(5) generator, and is a "momentum eigenstate" matrix satisfying [P^a, = 

Pda^(p) for complexified momentum Paa = AqAq. The corresponding currents for the plus 
helicity gravitons build up a self-dual background, and are given by: 

r+ = Tr([T,Q]DQ) (7.10) 

the minus helicity gravitons are specified by the related set of currents: 

r_ = Tr (^TQ^g) (7.11) 

which are acted upon by the 7+ generators. In an A/" = 4 supersymmetric theory these 
modes sit in different supermultiplets, and so can naturally be associated with different 
types of currents. Quite remarkably, the correlation functions of these currents reproduce 
the MHV graviton scattering amplitude of Einstein gravity in flat space-time. This cor- 
respondence is worked out in detail in [1], and forms the main piece of evidence that our 
large twistor matrix model gives rise to a space-time theory that contains gravity. 



7.3 Geometric Action 

In light of the detailed evidence that the defect sector of the matrix model contains the 
MHV sector of gravity, it is appropriate to ask whether the result can be extended to 
include non-MHV dynamics. Following the gauge theory lead, the logical place to look is in 
the hCS gauge sector and its coupling with the gravitational currents of the defect system. 
Indeed, we have already seen hints of a gravitational mode in the hCS system, in the form 
of an obstruction to decoupling the adjoint u{l) gauge transformationso 

To discuss the action of the adjoint u{l) in the large limit, it is helpful to refor- 
mulate the matrix multiplication of the matrix fields in terms of a Moyal product on the 
homogeneous coordinates: 

* = + (7. 12) 

for a (1, 1) bivector compatible with the symplectic structure induced by the Kahler 
form. The Moyal product for two functions / and g is, to leading order in g: 

f*g = f9 + Q''^ {djd^g - drfd^g) (7.13) 

where a and r are holomorphic and anti-holomorphic tangent bundle indices, respectively. 
Given a non-commutative U{1) gauge field Aa and a field Q in the adjoint, the covariant 

Though the physical context is different, for somewhat related discussion on non-commutative 
gauge theory and emergent gravity, see [53H57] and [58] for a review. 
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derivative is: 



VaQ = daQ + Aa* Q - Q * Aa. 



(7.14) 



In terms of the Moyal expansion, we have: 



(7.15) 



where we have introduced the modes: 



a'L = 2g''^drAa, 



(7.16) 



Note that what actually couples to the mode Q is the gradient of A rather than A itself. 
This again is an indication that the adjoint f/(l) behaves differently from other types of 
modes in the fuzzy system. 

The modes a and a transform as (0,1) forms valued respectively in T'^^'°^ and T(°'^\ 
i.e. the holomorphic and anti-holomorphic parts of the tangent bundle. A zero mode of Q 
is specified to zeroth order by the condition = 0. To first order, this is corrected to 
the condition d^Q — o^daQ = 0. By inspection, defines an almost complex structure, 
and the derivative Vq corresponds to the gauge field associated with complex structure 
deformations. This is of course in keeping with the expectation that twistor space is sensitive 
to complex structure deformations, and so a putative bulk gravitational theory will be of 
Kodaira-Spencer type. 

To further bolster this interpretation, we can also work out the (0, 2) component of the 
field strength J-"^^ = [Va, V^] . A lengthy but straightforward computation shows that the 
field strength takes the form: 



Here, {21, 21}^ denotes the Lie bracket for the (0, l)-form 21 taking values in the tangent 
bundle. In other words, the adjoint U{1) of the fuzzy system is naturally associated with 
the algebra of vector fields on twistor space. 

Given that the effective low energy theory describes physics derived from a holomorphic 
Chern-Simons field theory, it is natural to expect that the geometric action for the adjoint 
U{1) gauge field 21 can be packaged in the form of a BF type action that imposes the 
equation of motion J^*^"'^) = 0. So it is reasonable to write an action of the form 



where 21 is a (0, l)-form in the cotangent bundle which forms a pairing with the "color in- 
dices" of J-". The critical points of this action enforce the condition that the (0, 2) component 



where 21 = ^ {a^a + %i 





(7.18) 
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of the field strength is compatible with the complex structure of the geometry. Returning 
to our discussion between equations fl7.16p -( 17rT7|) . note further that to leading order in the 
gauge field fluctuations, only the almost complex structure deformation appears in a zero 
mode equation. Setting = 0, the field strength reduces to the Nijenhuis tensor for the 
almost complex structure. Vanishing of the field strength enforces the condition that it is 
integrable. This is of course quite familiar from earlier work on twistor string theory, and 
the appearance of a Kodaira-Spencer theory of gravity, as in [17]. BF type actions for the 
self-dual sector of supergravity have been considered in [5T1[59] (see also [60]). 

But as opposed to the twistor string, the coupling of the mode 21 to the defect system 
comes with additional structure, as reflected in the computation of MHV graviton ampli- 
tudes. Indeed, in the MHV amplitude, the glik^) charges which couple to G QS^'^\T'^^'^^) 
are invariant under longitudinal shifts Vaa ^aa + Vaa in the polarization tensor, which al- 
lows one to pick the gauge condition p ■ = 0. The positive helicity gravitons then describe 
special modes for which the (0, l)-form is Hamiltonian with respect to the infinity twistor 



7°^ [51]: 

a| = r'^a„/i^. (7.19) 

for h-^ a line bundle valued (0, l)-form. A similar redundancy condition holds for the dual 
mode 21 [51], and the corresponding physical mode h. These have the interpretation as the 
two helicities of a graviton. 

To give a bit more detail, recall that Einstein gravitons are given as elements h G 

i7"^(0(2)) and G i/^(0(— 6)), on PT', twistor space with a line at infinity removed. By 
contrast, the conformal graviton is instead represented by (0, l)-forms in the holomorphic 
tangent and cotangent bundle, i.e. C G i7^(T(^'°)) and C G i7'^(fi(^'°^). With no additional 
restrictions, C is simply the mode a of line f l7.16p . This is a linearized complex structure 
deformation for a Kodaira-Spencer theory of gravity on twistor space. To get back the 
usual modes of Einstein gravity, one introduces a redundancy for the C and C modes, so 

that only h and survive as physical excitations. As explained in [51], the precise 

link between the two is given by: 

where /"^ is the inverse bitwistor for Minkowski space. The redundancy is manifest for C, 
and for C is obtained through the conditions: 

C ^ C + dm + n (TTdn) (7.20) 
C^C + dx (7.21) 

where m is an element of fi(°'^)(0(-4)) and n is an element of f2*-°'^-'(0(— 6)). In addition. 
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X is an element of f2'^^'°^(0(— 4)), and the corresponding redundancy equation indicates C 
is represented by a cohomology class. 

On the 4D space-time, the exchange of the h/h pair between twistor lines (i.e. space- 
time points) then generates the usual gravitational potential for a spin two excitation with 
the strength of Newton's constant set by the ratio of an overall coefficient multiplying fl7.18p 
and the bulk defect coupling fl7.6p . It would be interesting to refine this discussion further 
and fix the precise value of Newton's constant. 

7.4 Einstein versus Conformal Gravity 

What is the low energy effective action for this putative theory of gravity? The leading 
order terms consistent with the symmetries of the system are: 

Sgrav = y d^Xy/^ (R - 2A) + (3 J d^Xy/^W^ + ... (7.22) 

where R is the scalar curvature, A is the cosmological constant, W is the Weyl tensor, and 
a and /3 are dimensionless couplings. In Einstein gravity, we would set a = l/(167rG'jvA). 
The case a = corresponds to conformal gravity. This latter theory contains a number of 
pathologies compared to ordinary Einstein gravity. For example, the kinetic term is quartic 
in the momentum, and the resulting spectrum of modes contains ghost-like excitations. 
This is also the gravitational sector of the twistor string theories of [TI IMIITT] . 

We will now list a set of arguments that indicate that in our case, rather than getting 
conformal gravity, the emergent gravity theory is described by the Einstein action with a 
Newton coupling set by a/GJv ~ ~ ^/ V^- 

• Graviton scattering: In principle, we can determine the effective action by computing 
all possible graviton scattering amplitudes which descend from the matrix model. So far, 
we have only computed the MHV amplitudes. Are these sufficient to distinguish the two 
theories? This is in fact a somewhat subtle question, because some features of Einstein 
gravity in (anti-)de Sitter space can be mimicked by conformal gravity with Neumann 
boundary conditions at spatial infinity [61] . 

MHV graviton scattering amplitudes describe the process of an incoming negative he- 
licity graviton bouncing off a self-dual space-time background. Self-dual backgrounds are 
solutions to both the Einstein and the conformal theory. So we can set up the calculation 
in both theories, by starting with a self-dual geometry that asymptotically looks like an 
(anti-)de Sitter space-time. Conformal gravity contains two propagating graviton modes, 
one that behaves like the ordinary graviton and a ghost mode. By picking appropriate 
asymptotic conditions, we can isolate the Einstein graviton. The conformal gravity calcu- 
lation then looks identical to the Einstein gravity calculation. The difference between the 
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two theories is that one theory has a dimensionful constant and the other a dimensionless 
one. The reason that we can compare the two directly is because we are considering the 
system on de Sitter space with finite curvature radius £^ ~ 1/A. 

For the sake of argument, let us suppose that the gravitational theory described by the 
matrix model is dominated by the conformal gravity term proportional to /3. The effective 
Newton's constant would then fixed by the coefficient (3 and the de Sitter radius ^iT\ : 

If we were to match our space-time theory to this prescription, we would need to set 
A ~ l/i"^ and Gn ~ ~ ^^/N which would imply that P N. In the large iV/ffat space 
scaling limit, this looks absurd. Indeed, conformal gravity (with any finite couphng /3) does 
not have a non-zero MHV graviton scattering amplitude in fiat spaced Further, explicit 
calculations of graviton scattering in twistor string theory vanish for Einstein gravitons [62] . 
Our double scaled matrix model does have finite fiat space scattering amplitudes. 

• Absence of conformal symmetry at finite ipf. Conformal symmetry is explicitly broken 
in the double scaled matrix model, both due to the presence of a UV cut-off ipi and due 
to the introduction of the infinity bi-twistor. In the non-commutative twistor theory, only 
the compact SU{4:) subgroup of the complexified conformal group S'L(4, C) is unitarily 
realized. Moreover, the defect action, which involves the infinity twistor, further breaks the 
symmetry to 5*0(5). The 5*0 (5) generators therefore play a distinguished role among the 
currents of the defect system. The (MHV manifestation) of the graviton modes are made up 
from these 50(5) currents - specifically, those associated with the hermitian 'translation' 
generators (17. 8p . The hermitian generators (17. 8p are precisely the ones that are compatible 
with the boundary conditions [HI] that selects the Einstein gravitons from the conformal 
gravity fiuctuations@ 

• CSW correspondence: As explained in section 6.1 and in [1], integrating out the 
defect modes produces an effective action that, in the continuum limit, matches with the 
MHV effective action in [151 SS], which (by design in [13 US], but derived here) evades 
the troubling appearance of conformal gravity amplitudes, that plague the twistor string 
proposal of [3117] . The appearance of gravity for us indeed has a different geometric origin. 

Taking /3 ~ — )• oo amounts to taking a classical limit in conformal gravity. One could sus- 
pect that this classical limit amounts to a decoupling limit, in which the ghost mode of conformal 
gravity can be consistently decoupled from the rest of the theory. If this were possible, the left 
over graviton mode would propagate and interact like an ordinary graviton in Einstein gravity, 
with a finite Planck length equal to £pi. 

^'^Note also that, by starting from the with finite radius £, we are able to define hermitian 
translation generators ()7.8p . and avoid the non-unitarity that plagued the twistor string theory 

of HZ]. 
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in which the Einstein gravitons immediately play a distinguished role. 

Suppose the theory does indeed give rise to Einstein gravity, as we have argued. What 
kind of theory is it, and does the matrix model - or rather its UV parent theory, the 
holomorphic Chern-Simons theory with flux - provide a UV complete description? With 
regard to the flrst question, indications are that the IR theory will take the form of A/" = 4 
supergravity, possibly coupled to A/" = 4 gauge theory. We have not found a clear condition 
that limits the possible rank of the gauge group. This suggests that the answer to the second 
question should be: No, the matrix model (and even its embedding in the topological B- 
model string theory) is not UV complete, but should be viewed as the low energy effective 
description of a more complete theory. To illustrate how such an embedding may arise, we 
now describe a potential realization of the matrix model in the physical superstring. 

8 Embedding in Superstring Theory 

In this section we propose an embedding of the twistor matrix model in a brane construction 
in superstring theory. In particular, this will establish a concrete UV completion for our 
system. The basic setup we consider is given by a DO-brane bound to a stack of eight-branes 
in type IIA superstring theory. We work in ten non-compact space-time dimensions, which 
we write as M^ume x x ]R_|_. The eight-branes are assumed to be parallel, fllling M^ume x 
and sitting at various points of M^. The DO-brane we shall be considering will be bound to 
one stack of eight-branes, but will be free to move in the worldvolume M^jme x M^. The basic 
idea is that when a suitable flux is switched on along the worldvolume of the eight-branes, 
we obtain four copies of the supersymmetric harmonic oscillator. The statistical mechanics 
of this system is studied by compactifying on a thermal circle of radius /3. We argue that 
the many body problem for a gas of 0-8 strings at low temperature realizes the twistor 
matrix model. 

8.1 Brane Construction 

Our main interest will be in the worldvolume theory of a single DO-brane probing a stack 
of Nds D8-branes coincident with an OS-plane. However, to frame our discussion, let us 
flrst briefly review the worldvolume theory of Ndq DO-branes, flrst in the absence of other 
branes, and then in the presence of eight-branes. 

Consider flrst the case of Ndq DO-branes in ten flat space-time directions. This is a 
system which preserves sixteen real supercharges, and is described at low velocity by a 
quantum mechanics with (8, 8) supersymmetry. For N^q DO-branes, we have a UINdq) 
gauge theory with gauge fleld Ado- All of the 0-0 strings transform in the adjoint of 
U{Ndo)- The bosonic matter content consists of eight real scalars for 1 = 1, .., 8 which 
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describe motion parallel to the D8-brane and an additional real scalar X± describing motion 
transverse to the D8-brane. In addition there are sixteen fermions, which transform as a 
Majorana-Weyl spinor of so(9, 1). Under the decomposition so(9, 1) D so{8) x so(l,l), 
these fermions transform in the 8^ and 8c of so{8). We denote these two sets by 5*^ © Sa' 
for A, A' = 1,...,8. 

Introducing a stack of D8-branes leads to additional dynamics for the system. Due to 
their high codimension, adding D8-branes introduces a number of subtleties. For example, 
at finite distance from the D8-brane, the dilaton blows up. Consistent treatment of the 
system then requires introducing 0-planes to prevent this. When the tadpole is not locally 
cancelled, there will also be a Romans mass on one side of the D8-brane |63j . 

The effective theory of the DO-brane in the presence of Nds D8-branes has been treated 
previously in the literature, see for example [01]. This system preserves a chiral (8,0) 
super symmetry. In addition to the 0-0 strings, the worldvolume theory now contains an 
additional 0-8 string charged in the {N^s, Ndq), which is a complex fermion we denote by 
X- This mode is a singlet under the (on-shell) supersymmetry transformations, which is 
possible due to the low dimension and chiral nature of the supersymmetry. The resulting 
D0-D8 system is BPS. In particular, the DO-brane experiences no force in the direction 
transverse to the D8-brane. See |6l] for further discussion. 

The system we shall mainly be interested in for our purposes is the case of N^s D8- 
branes coincident with an 08-plane. Local tadpole cancellation then fixes N^s = 16, that 
is, the maximal gauge group is 5*0(16). Let us note that £^8 can also be realized, though 
non-perturbatively. The presence of the orient if old plane affects the 0-0 and 0-8 strings as 
well. The DO-brane gauge group is now ©(A'^dq)!^ The 0-0 strings, which were initially 
in the adjoint of U{Noo) are now two index representations of 0{N£,q). The modes in the 
two index symmetric representation descend from © Sa, while the modes in the two 
index anti-symmetric descend from A^q (B Sa' © X±. The 0-8 strings now transform in 
the {Noo, Nds), which is a real representation. To maintain a holomorphy convention, we 
shall often decompose SO^Nds) ^ U{N£,%I'2). The 0-8 string content then consists of a 
vector-like pair of complex fermions x®X, which are subject to the reality condition = X- 

Consider now the special case of a single DO-brane. When N^o = 1, we project out the 
mode describing motion transverse to the D8 stack, and also eliminate half of the fermionic 
modes. The gauge group is also in this case reduced to 0(1) = The effective Lagrangian 
is (in units where the DO-brane mass is unity): 

Le// = ^dtX'dtX' + '-SAdtSA + ixDx + xDx) (8-1) 

^^As explained for example in [65] the reason the gauge group is 0{Ndq) rather than SO{N£)q) 
is because upon compactifying on the transverse circle 5*^, there are Z2 valued Wilson line con- 
figurations available in the T-dual description. 
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where the covariant derivative acts on the fermionic modes as: 



Dx={dt + igAr>^)x (8-2) 
Dx={dt-igAm)x (8.3) 

where Aos is the puUback of the D8-brane gauge field to the DO-brane worldvolume, with 
g the gauge couphng of the D8-brane theory. In the above, we have speciahzed to the case 
of a gauge field activated in the U{Ndq/2) subgroup. 

This system describes a DO-brane which is free to move inside of the D8-brane. In 
the quantum theory, the conjugate momenta to are 11^ = dtX^ . This results in the 
commutation relations: 

[X', U'] = {Sa, Sb} = 6ab, {xk, xj} = hi (8-4) 

where k (resp. I) denotes an index in the fundamental (resp. anti-fundamental) of U{N£,s/2). 
Observe that there is a degeneracy of ground states for the system. These correspond to 
the possible positions of the DO-brane in the superspace R®'^. We will soon switch on a fiux 
which identifies a canonical choice of complex structure, in which case we have C^'^. 

To break this degeneracy in the ground state energies we now switch on a magnetic fiux 
which takes values in the overall U{1) C U{Nj:i^/2). The effective energy scale for the flux 
is a tunable parameter of the D8-brane worldvolume theory. Since the DO-brane is more 
accurately thought of as a superparticle, it is appropriate to consider a background gauge 
fleld configuration: 

Ads = Fi'^X'dtX' + F^iS^S"" (8.5) 

where here, F^""^ and F^^^ refiect the presentation of the fiux in the 8^, and 8^ of so(8). By 
triality, we can take an embedding of su{4) x m(1) c so{8) so that the 8^, and 8^ decompose 
as: 

so(8) D sm(4) X u{l) (8.6) 
8, ^4+i + 4_i (8.7) 
8,^4^1 + 4+1 (8.8) 

The original basis of modes can now be traded for the usual coordinates of supertwistor 
space, © In terms of this basis, the puUback of the gauge field is: 

Ads = F^-p [z^^d^W - Z^dtZ^^ + F.-j^j^ (8-9) 
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The profile of tliis flux is taken to be: 



F = i^dz^ A dz^ + i^dz^ A dz^ + i^dz^ A dz^ + i^dz^ A dz\ (8. 10) 

2 2 2 2 ^ ^ 

In general, this choice of flux will break supersymmetry on the D8-brane, though the DO- 
brane may still experience an effective supersymmetric quantum mechanics!^ 

Let us now study the quantum mechanics of the endpoint for the 0-8 string which is 
attached to the D8-brane. Evaluating in the background where the 0-8 string is present, 
the Lagrangian is: 

L,ff = ^dtX'dtX' + '-SAdtSA + igFl'fx'dtX' + igF^is^S''. (8.14) 

Quantization is now of the standard type for a superparticle moving in a background mag- 
netic flux. This is described by four supersymmetric harmonic oscillators, which by abuse 
of notation we label as and ip^. These modes have the same supercommutators as fuzzy 
supertwistor space, with characteristic frequencies: 

Ua = (8.15) 

Here we have rescaled the coordinates to include the explicit dependence on the DO-brane 



^^Following |66H68| . to study the unbroken super symmetries in the presence of this flux, intro- 
duce the block diagonal matrix M in the 8^ representation (see e.g. [6711691 ^ 



4 

5 
=1 



M= (BMi, Mi 



cos 2TTVi sin lirvi 
— sin 27rvi cos 27rvi 



.11) 



where each Vj is given by: 



exp 2TTivj = ^^^^^ . (8.12) 
1 — ijj 



The condition for unbroken supersymmetry is that in the 8c representation, p{—M) has a unit 
eigenvalue. The amount of supersymmetry which is preserved is dictated by the maximal subalge- 
bra of it(4) C so{8) which commutes with p{—M). A necessary condition for some supersymmetry 
to be preserved is: 

vi + ... + V4eZ. (8.13) 

For a generic choice which retains some supersymmetry, this leaves a system with only two real 
supercharges. The case we shall eventually specialize to is given by taking all fj = f- In this case, 
we see that generically, all supersymmetry will be broken. However, in the limit / — )■ +oo, we see 
that all eight supercharges are retained. This corresponds to the limit Vi — )■ 1/2. 
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Figure 3: Depiction of the DO-brane bound to a stack of N,. D8-branes and an 08-plane in 
the presence of a magnetic flux threading the D8-brane. The 0-8 strings are attached to the 
DO-brane, and the other end can wander inside the D8-brane. Switching on a background 
flux through the D8-branes, a gas of 0-8 strings descend into Landau levels graded by 
integers M > 0. Here, we have depicted a low temperature configuration which realizes 
fuzzy supertwistor space at level = Mmax- 

mass. The effective Hamiltonian is: 

4 4 

Heff = Y^^aZiZ^ + Y^U^^i^' (8.16) 

a=l i=l 

where the zero point energy cancels out. In what follows we focus on the special case of 
iu = oja = for all a and i. Indeed, in this case is proportional to the Hamiltonian 
constraint operator Hq of equation fl3.18p . 

An important feature of this system is that it admits a decoupling limit, where we take 
/ to be large in string units, while keeping the ratio f /moo small in string units. This 
is basically the same type of limit studied in [72] for non-commutative Yang-Mills theory. 
Here, this zero slope limit ensures that the low energy dynamics are well-described by this 
effective Hamiltonian. 

The modes at level M are described by the Hilbert space of fuzzy points for CP^'^, 
?^j.p3|4(M), so there is a degeneracy of Km at each M. Thus, the phase space available 
to a 0-8 string is now a direct sum of fuzzy supertwistor spaces. Finally, at each level, 
there is an additional overall degeneracy, because the fermionic 0-8 string transforms in 
the fundamental of U{NDg/2). We denote the effective Nc = Nds/2. See figure E] for a 
depiction of this system. 
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8.2 Statistical Mechanics 



Let us now study the statistical mechanics of a large number of 0-8 strings. We compactify 
the temporal direction on a circle of circumference /3 = 1/T and take thermal boundary 
conditions for all modes. Though there are some subtleties in such stringy systems at high 
temperature, i.e. small thermal radius, we shall mainly confine our discussion to the low 
temperature limit. 

In the grand canonical ensemble, the chemical potential for the 0-8 strings is controlled 
by a flat connection from the D8-brane gauge field. The full gauge field configuration on 
the D8-brane is then given by the magnetic flux considered earlier, and a holonomy from 
the component which points along the thermal circle direction. The energetics of the 0-8 
strings are unaffected by including the contribution from the flat connection, though it 
does introduce a chemical potential for the system. Integrating over the thermal circle, we 
obtain the fugacity ( for the grand canonical ensemble: 




C = exp / Ans\. (8.17) 



The chemical potential is then given by /i = TlogC- 

Consider first the zero temperature limit, with a fixed number of 0-8 strings which we 
denote by Ntot. The phase space for each particle is a non-commutative C^'^. Since the 0-8 
strings are fermionic, they will start to fill up a Fermi surface. We will mainly be interested 
in the case where the Fermi surface forms a complete shell up to a maximal level Mmax = 
This corresponds to taking a fixed number of states: 

ax 



A/=0 

Each energy shell corresponds to the Hilbert space of a fuzzy CP'^'^ at level M. As we 
increase Mmax, we approach a large limit. The states fill out entries in the direct sum: 

ax , , , , 

niNtot) = H^(m {N) ®N,= © 'H„-,i4{M) ® N, (8.19) 

M=0 

where the tensor product with the factor of A^^^ reflects the fact that each particle transforms 
in the fundamental of U{Nc). Up to this factor, we can identify this Hilbert space with the 
fuzzy four-ball encountered previously in section HJ 

The length scale £ is set by the inverse characteristic frequency of the system u: 

ir^uj-\ (8.20) 
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We can also introduce an anisotropic scaling for the D8-brane magnetic flux, and thereby 
include the parameter 7 of the infinity bitwistor. 

Returning to the case of uniform energy level spacing in the phase space, we would now 
like to model the effective dynamics of the many body system at finite temperature. It 
is helpful to work in terms of fields, which describe the collective excitations. There are 
two types of excitations which we can identify, corresponding to motion within the Fermi 
surface, and normal to it. 

Excitations above the Fermi surface are 0-8 strings, which are characterized by a field 
Q{Z\ip\ Z,ip) in the fundamental of U{Nc). Holes are 8-0 strings, and are characterized 
by a field Q{Z^ , ip\ Z, ip) in the anti-fundamental. Moving the location of the DO-brane, we 
see that there are also four complex scalars. These can be viewed as bound states of the 
8 — © — 8 strings. If the particles were free to move inside of the D8-brane, we would 
have four such complex scalars ^"(Z"'", ■?/;''", Z, These modes define the transformation 
of the states within the Fermi surface, and are therefore maps ^{Ntot) — ^ 'H{Ntot + !)• In 
other words, they are fuzzy (0, l)-forms acting on 'H^{4\4){N) which are also in the adjoint 
representation of U{Nc). Since we are interested in the low energy dynamics confined to 
near the Fermi surface, there is an additional constraint: 



so there are only three independent bosonic degrees of freedom. 

We fix the effective action via symmetry considerations. The kinetic operator for the 
particle/hole pairs corresponds to an excitation in the direction normal to the Fermi surface. 
In a local patch where we write CP^'^ ^ C^l"^ x CP^I°, this is in the direction of the bosonic 
twistor line. From our earlier analysis of the defect kinetic operator, we know that a 
derivative operator which is local on the "space-time" directions of the twistor fibration is 
the kinetic operator V. Symmetry consideration then dictate the kinetic operator for the 
particle/hole pair at low temperatures: 



where the form of the derivative operators is fixed by gauge invariance. Here, the defect 
action is integrated over the entire four-ball of the phase space. Excitations near the top 
of the Fermi surface cost less energy, and so at low energies, we get back the integral over 
just the top layer, given by the CP^'^ at level = Mmax- 

For the bulk modes, there is a single topological term we can write down which is 

Although it is tempting to also include four complex fermionic modes {Z^ , Z , ip) , the 
geometric interpretation of these modes is less clear, as we cannot "give a vev" to a fermionic 
direction. 



ZlA'^ = 0. 



(8.21) 




(8.22) 
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consistent with the symmetries of the system; it is the holomorphic volume of the phase 
space in the supermanifold C^'^ occupied by the 0-8 strings: 

^buik = Ci A 1th(n,.^) [s^p.sD'-dPd^D')^, (8.23) 

where = + A"" and Ci is the RR one-form potential. The coupling constant of the 
matrix model is obtained by integrating over the Euclidean thermal circle. Taking Ci to 
be a constant one-form and integrating over the circle yields an overall coupling. In other 
words, the action becomes: 

5buik m = ^Tr^(iv..o [s^p.sD'^DPD^D')^, (8.24) 

Hence, the low temperature dynamics of a large number of 0-8 strings reproduces the basic 
form of the twistor matrix model. 

By embedding the matrix model in the physical superstring, we have equipped the 
model with a proposed UV completion. Consistency of the UV completion leads to some 
interesting restrictions on the form of the low energy effective theory. For example, it 
imposes a sharp upper bound on the rank of the gauge group. Note, however, that this 
restriction is milder than what is required in conformal A/" = 4 supergravity for anomaly 
cancellation (see [73] for a review). A further interesting feature of this construction is that 
it is formulated with ten non-compact directions. Indeed, the embedding of the twistor 
matrix model into the superstring does not involve a traditional compactification at all0 

Using this perspective, we can ask how various features of the 4D space-time are encoded 
in the gas of 0-8 strings. First note that the double scaling limit for realizing a fiat space 
theory is quite natural in this setup; we can scale N — )■ oo, but simultaneously decrease the 
energy spacing between the Landau levels. The currents of the matrix model are bilinears 
in the QQ modes, corresponding to the creation of particle/hole pairs. If we are at low 
energies, all of these states will be created near the top of the Fermi surface, and at the same 
Landau level. However, as we increase the available energy, holes can be created further 
down in the Fermi sea, and can also jump between neighboring Landau levels. Thus, a 
continuous extra dimension opens up in the direction normal to the Fermi surface. It is 
tempting to identify this with a holographic RG scale. At even higher excitation energies, 
the particles are no longer close to the Fermi surface, and simply move throughout the 
fuzzy C^l^ phase space. Note that even in this case, however, there is an upper bound on 

^^In string theory the usual situation is that without a compactification of the "extra dimen- 
sions" one obtains a 4D gauge theory which is decoupled from gravity. For example, a D3-brane 
probing a non-compact Ricci flat manifold will achieve this. Here, the interpretation of the ten 
dimensions for the 4D theory is different, and gravity emerges along with the space-time. 
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the available excitation energies. 

Though at finite the continuum theory interpretation is less evident, it is also the 
case where we have a truncated Hilbert space for Euclidean de Sitter. In this case, the 
jump in Landau levels is discontinuous. A further curious feature is that the effective ratio 
^/ ^pi ~ is modified at high energies so that states would appear to experience distinct 
values of the cosmological radius. One can also consider the high temperature regime ~ 
i.e. high compared to the level spacing - of the gas of 0-8 strings, but below the Hagedorn 
temperature. This is still a well-defined system, but now the Fermi surface itself will begin 
to break apart due to thermal fiuctuations. Clearly, it would be interesting to study this 
set of issues further. 

8.3 Dual Descriptions 

The D0-D8 bound state we have been considering has a number of dual formulations, which 
provide additional insight into the dynamics of the system. Compactifying on a circle S\_ 
in the direction transverse to the eight-branes, we obtain a T-dual description in type I 
string theory, where the DO-brane is now a Dl-brane. This is in turn dual to an Fl-string 
of IIB via an S-duality. Using type I/heterotic duality, we can map this to a heterotic string 
wrapping Mtj^e x S\. Observe that the projection on the mode content we have observed 
is quite natural on the heterotic string side. There is also an intriguing connection to the 
M = 2 string. Indeed, as advocated for example in [71], the worldvolume theory of the 
Dl-brane is naturally viewed as a four-dimensional theory on M^'^ which has been reduced 
along a null vector. This suggests a more direct connection between M = 2 strings and the 
physics of supertwistor space which would be interesting to develop further. 

We can also use these dual formulations to realize more intricate gauge theories. The 
basic idea is to consider the low energy quantum mechanics of a DO-brane, now bound to 
a more general configuration of eight-branes and lower dimensional branes. Compactifying 
on a circle, this type of configuration is T-dual to a configuration of intersecting seven- 
branes in IIB. Going to strong coupling, we can study the resulting dynamics in terms of 
an effective heterotic string probing an appropriate flux background. 

As an example, consider three stacks of branes which we denote as S, B' and B" which 
wrap a common Mfj^e x I^^- We assume B describes a D8-brane, as in our previous discus- 
sion. Though it is beyond the scope of our present discussion to provide a full description 
of this system, we can sketch how additional matter sectors could arise. In addition to the 
— i3 strings, there will now be — i3 and — B" strings. Excitations above the Fermi 
surface can therefore be of different particle/hole types: They can correspond to pairs of 
S — © — S modes, which we previously identifled with the modes of an A/" = 4 supermul- 
tiplet in the adjoint of U{Nc). In addition, we see that there are new effective modes, such 
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as ;B — © — i3' which are in the fundamental of the emergent U (Nc). Observe that gauge 
invariance also dictates the form of the possible interaction terms of these composites. 

Note also that the effective dynamics is not limited to unitary gauge groups. Indeed, 
using the dual description in the heterotic string, there is an effective worldsheet description 
available in the strong coupling limit of the IIB configurations, so we can also consider E- 
type gauge theories. This suggests that phenomenologically relevant gauge theories may 
also be engineered in this type of setup (by embedding in a GUT group). It would be of 
interest to provide a more explicit constructions along these lines. 

9 Conclusions 

In this paper we have provided evidence that A/" = 4 SYM theory on an is dual to a 
large N matrix model defined over non-commutative twistor space. The matrix model is 
the theory of the lowest Landau level for holomorphic Chern-Simons on CP'^'^ expanded 
around the Yang monopole configuration. The matrix model incorporates the symmetries 
of the original space-time theory, and moreover, contains a natural class of spin 1 and spin 
2 symmetry currents. Quite remarkably, in the fiat space limit of the matrix model, the 
correlators of these currents correctly reproduce MHV gluon and graviton scattering pQ. 
This is a highly non-trivial feature, and provides evidence that the low energy limit is 
described by Einstein gravity. We have also presented a proposed UV completion of the 
model based on the low energy dynamics of a D0-D8 bound state. In the remainder of this 
section we discuss various open directions. 

At a pragmatic level, it is important to check that the match to the 4D theory continues 
to work at loop level in the gauge theory. Indeed, loop level contributions from conformal 
graviton states is a significant hurdle for the general twistor string program. This would 
also provide another probe into the nature of the gravitational theory which is coupled to 
the gauge theory. 

The most intriguing feature of the twistor matrix model is that it seems to contain an 
emergent gravitational sector. Though there are strong constraints on emergent theories of 
gravity, the basic setup of the matrix model violates some assumptions of the Weinberg- 
Witten theorem [75] . For example, the graviton and the space-time emerge simultaneously 
in a double scaling limit. Furthermore, the theory is defined as the infinite radius limit of a 
theory with positive curvature. Additionally, it is far from evident that our model possesses 
a local gauge invariant stress energy tensoro Note, however, that although there is a short 

^^One could ask whether the model comes with additional higher spin currents, which did appear 
in [9]. Here, note that the interactions involving our spin 2 excitation are suppressed by the small 
parameter ipi, and so in the large hmit, it is natural to expect any higher spin excitations to 
effectively decouple. Even in the finite setting, these excitations are more appropriately viewed 
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distance cutoff the low energy theory will not contain Lorentz violating operators. This 
is simply because the model respects so(5) and its Wigner-Inonii contraction in the ffat 
space limit. To further test our conjecture, it would be interesting to see whether other 
features of quantum gravity such as black hole formation, or a more direct account of 
holographic entropy bounds are also realized. 

We have also seen that fuzzy twistor geometry and the matrix model may naturally arise 
from a D0-D8 bound state. It is likely that this construction generalizes to situations with 
less supersymmetry and more realistic matter content. Observe that all ten of the string 
space-time dimensions are necessary to realize the 4D dual theory. Thus, four space-time 
dimensions is a rather special part of the construction. 

Finally, as the 4D theory is formulated on a four-sphere, it is natural to analytically 
continue the theory to de Sitter space. In particular, the large but finite version of the 
matrix model suggests an intrinsic link between the de Sitter radius and the Planck length 
via ~ N. In other words, the cosmological constant is always parametrically small 

in Planck units. The finite number of degrees of freedom at finite is also suggestive of a 
holographic theory of de Sitter space [80] . 
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